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,  ■  ABSTRACT 

work  we  stu^jr  subharmonic  solutions  near  an  equilibrium  point  for  -the 
Hamiltonian  system  , 

- - z  =■  J A{t)z JHj{z,t) 

where  A(t)  is  a  matrix,  Hz{z^t)  =  o(|  z  ()  and  both  A  and  H  are  periodic  in  t. 

On  the  linear  part  of  the  system  we  impose  a  condition  expressed  in  terms  of  its 
symplectic  invariants.  The  higher  order  term  is  assumed  to  be  superquadratic  near 
the  equilibrium  point,  and  we  show  that  this  condition  can  be  reduced  to  the  center 
manifold. 

We  transform  the  Hamiltonian  system  to  a  variational  problem  and  we  apply  a 
minimax  argument  to  find  critical  points.  \  ,  .i  -  - 
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SUBHARMONIC  SOLUTIONS  NEAR  AN  EQUILIBRIUM  POINT 
FOR  HAMILTONIAN  SYSTEMS. 


§0  INTRODUCTION. 

In  this  work  we  study  the  existence  of  subharmonic  solutions  near  an  equilibrium 
point  for  Hamiltonian  systems 

z  =  JHz(z,t)  (3-1) 

where  z  =  (p,  5)  G  R^”,  J  denotes  the  standard  symplectic  structure  in  and 
denotes  the  p^lrti^Ll  derivative  of  H  with  respect  to  z. 

We  eissume  that  the  Hamiltonian  H  satisfies: 

(HO)  H  is  near  z  =  Q, 

(HI)  /r(0,<)  =  0,  /f^(0,<)  =  0  V<€R,  and 

(H2)  H  is  T-periodic  in  the  t-variable. 

Hypothesis  (HO)  and  (Hi)  allows  us  to  write  the  Hamiltonian  as 

H{z,  t)  =  ^{A{t)z,  z)  +  H{z,  t)  (0.2) 

where  A{t)  denotes  the  Hessian  matrix  of  at  z  =  0  and  H{z,t)  =  o(|zp)  represents 
the  higher  order  terms  of  the  Hamiltonian. 

We  then  can  rewrite  system  (0.1)  as 

z  =  JA{t)z  +  JHz{z,t),  (HS) 


and  z  =  0  represents  an  equilibrium  point  of  (HS). 

The  problem  of  finding  subharmonic  solutions  near  z  =  0  for  (HS)  has  been  the 
object  of  a  considerable  amount  of  work  (See  [2],  [3],  [7],  [11],  [12],  and  references 
therein).  The  first  results  are  due  to  Birkhoff  and  Lewis  in  1933,  [3].  In  1968  Harris  [7] 
gave  a  more  general  version  of  the  results  of  Birkhoff  and  Lewis.  And  in  1977  Moser 
[11]  presented  an  improved  version  of  those  mentioned  above. 

The  work  in  [3]  emd  [7]  requires  a  strong  nonresonance  condition  on  the  Floquet 
exponents  of  the  linearized  system.  This  allows  reduction  of  the  Hamiltonian  to  the 
Birkhoff  normal  form.  With  an  additional  eissumption  of  nondegeneracy  on  the  higher 
order  term  they  prove  the  existence  of  subhzu’monics.  r 

Moser  in  [11]  restricted  the  nonresonance  condition  only  to  the  Floquet  exponents 
on  the  imaginary  axis  by  reducing  the  problem  to  the  center  manifold.  The  nondegener¬ 
acy  condition  on  the  higher  order  terms  was  also  imposed  only  on  the  center  manifold. 
Moser  presented  his  results  working  in  the  context  of  symplectic  diffeomorphisms  rather 
than  directly  with  (HS).  The  differentiability  hypothesis  on  the  map  was  also  reduced 
to  the  class  C^.  Piatributiog/ 
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The  work  of  Rabinowitz  in  [12]  gave  another  direction  in  the  methodology  to  attack 
the  problem.  By  using  minimax  techniques  of  a  global  nature  he  proved  the  existence 
of  subharmonics.  He  assumed  that  the  linearized  system  has  all  its  Floquet  exponents 
on  the  imaginary  axis,  and  different  from  each  other,  and  that  the  higher  order  term 
is  superquadratic  in  an  appropriate  way.  He  also  obtain  the  result  when  the  linearized 
system  is  trivial,  i.e.  A(t)  =  0. 

Recently,  Benci  and  Fortunato  [2]  gave  weaker  conditions  on  the  Floquet  exponents 
of  the  linearized  system,  keeping  the  higher  order  term  H  superquadratic. 

Comparing  the  results  of  [3],  [7],  and  [11]  with  those  in  [2]  and  [12],  we  see  that  in 
the  latter  the  condition  on  the  linear  part  is  relaxed  by  assuming  a  stronger  condition 
on  the  higher  order  term. 

In  this  work  we  generalize  results  in  [2]  and  [12].  In  the  first  place  we  relax  the 
conditions  on  the  Floquet  exponents.  In  the  second  place  we  reduce  the  system  to  the 
center  manifold  and  we  require  the  higher  order  terms  of  the  reduced  Hamiltonian  to 
be  superquadratic. 

A  more  precise  description  of  the  results  follows.  Let  us  consider  the  Lneaxized 
system 

i  =  JA{t)z.  (0.3) 

After  a  coordinate  transformation  (0.3)  can  be  written  as 


?  =  JBz  (0.4) 

where  H  is  a  constant  symmetric  matrix.  The  matrix  JB  can  be  further  transformed 
into  a  canonical  form  where  only  symplectic  cheinges  of  coordinates  are  allowed. 

It  is  the  nature  of  the  canonical  fcrm  oi  JB  that  will  determine  the  possibility  of 
solving  (HS)  for  subhaxmonics. 

We  do  not  expect  to  find  subharmonic  solutions  in  the  case  all  Floquet  exponents 
are  outside  the  imaginary  axis.  See  examples  in  [2]  and  [11].  Thus  we  will  assume  that 
there  is  at  least  one  purely  imaginary  Floquet  exponent.  However  ,  as  we  will  see  in 
Section  3,  this  is  not  sufficient. 

In  Section  2  we  will  define  a  Floquet  exponent  of  positive,  negative  or  indefinite 
type.  This  is  a  symplectic  invzu-iant  related  to  a  perturbed  eigenvalue  problem  for  the 
blocks  of  the  canonical  form  of  JB. 

In  terms  of  this  notion,  we  assume 

(H3)  The  linear  system  (0.3)  has  at  least  one  purely  imaginary  Floquet  exponent  that 
is  not  of  negative  type. 

For  the  higher  order  term  H  we  will  assume 
(H4)  There  are  constants  r>0,  c>0, 5  +  2>p>g>2  that 

iHA^.t),z)>qHiz,t)>c\z\P, 
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VzeR^”,|z|<r  and  VfeR- 
Then  the  following  theorem  holds. 

Theorem  I. 

Suppose  (HO),  (Hi),  (H2),  (H3)  and  (H4)  hold.  Then  there  is  an  integer  N  so 
that  for  every  k  >  N  there  exists  a  nontrivial  solution  zjt  of  (HS)  with  period  kT,  and 
Mr=  ^  converges  to  zero  as  k  oo  uniformly  in  C^(R,  R^”). 

Remark  0.1.  Benci  and  Fortunato  in  [2]  considered  the  following  hypothesis  on  the 
linezirized  system  (0.3). 

(H3’)  Every  symplectic  matrix  sufficiently  close  to  X{T)  heis  at  least  one  eigenvalue  on 
the  unit  circle. 

Here  X(t)  represents  the  fundamental  matrix  of  (0.3),  i.e.  X(t)  satisfies 


X(0)  =  hn 

(0.5) 

X{t)  =  JA(t)X(t). 

(0.6) 

We  will  show  that  (H3’)  implies  (H3),  but  the  converse  is  not  true.  See  Section  3.3. 

Remark  0.2.  In  [12]  results  similar  to  Theorem  I  were  given  in  the  case  i4(^)  =  B  is 
constant  and  either  B  =  0  or  JB  has  purely  imaginary  eigenvalues,  all  different.  We 
will  see  that  both  cases  satisfy  hypothesis  (H3).  We  note  that  B  =  0  does  not  satisfy 
hypothesis  (H3’)  of  Benci  and  Fortunato. 

Remark  0.3.  In  the  election  of  the  names  Floquet  exponent  of  positive,  negative  or 
indefinite  type  we  want  to  reflect  the  possibility  of  finding  subharmonics  for  (HS)  with 
the  higher  order  term  H  satisfying  (H4)  or  the  opposite; 

(H4-)  There  are  constants  r>0,  c>0,g  +  2>p>g>2  such  that 

{H,(,z,t),z)  <  -qH(z,t)  <  -c\z\P 

Vz6R^",|zl<r,VfGR, 

or  both.  See  Remark  3.4. 

Remark  0.4  Our  hypothesis  (H4)  is  slightly  more  general  that  the  corresponding  one 
in  [2]  and  [12].  They  assume  p  =  q.  This  will  be  needed  in  the  results  obtained  on  the 
center  manifold.  See  Section  4. 

Remark  0.5.  If  we  further  assume  that  0  is  not  a  Floquet  exponent  for  (0.3)  then 
using  an  argument  given  in  [2],  we  can  prove  that  the  solutions  given  by  Theorem  I 
have  minimal  period  kT,  for  k  large  and  prime.  This  hypothesis  requires  that  the 
only  T-periodic  solution  of  (0.3)  be  the  trivied  one.  It  will  be  violated  if  (HS)  is  the 
variational  equation  about  a  periodic  solution  of  an  autonomous  system. 
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Theorem  I  gives  a  result  generalizing  earlier  works  by  relaxing  conditions  on  the 
linearized  system.  The  second  goal  of  this  work  is  to  generalize  condition  (H4),  that 
has  to  do  with  the  nonlinear  term  H. 

By  assuming  extra  differentiability 

(H5)  H  is  near  z  =  0, 

we  reduce  system  (HS)  to  the  center  manifold.  Specifically  v/e  will  obtain  a  Hamiltonian 
system 

w  =  JB{t)w  +  (0.7) 

where  the  linearized  system 

w  =  JB{t)w  (0-8) 

has  all  Floquet  exponents  on  the  imaginary  axis. 

We  assume 

(H6)  There  are  constants  c>0,  r>0, 5  +  2>p>9>2  such  that 

(Ku>{w,t),w)  >  qK{w,t)  >  c|i(;|^ 

Vuj€  <r,V<6  R. 

where  2m  is  the  dimension  of  the  center  manifold  of  (HS). 

An  application  of  Theorem  I  will  give 

Theorem  II. 

If  (H5),  (Hi ),  (H2),  (H3)  and  (H6)  hold  then  there  exists  an  integer  N  so  that  for 
every  k  >  N  there  exists  a  nontrivial  solution  Zk  of  (HS)  with  period  kT  and 
converges  to  zero  as  k  oo  uniformly  in  C^(R,R^"). 

Remark  0.6  Since  hypothesis  (H6)  is  in  general  hard  to  check  we  are  interested  in 
giving  conditions  on  the  original  function  H  that  imply  (H6). 

Roughly  speciking,  if  the  superquadratic  condition  in  (H4)  is  satisfied  for  every  z 
on  the  center  manifold  then  (H6)  holds.  A  precise  formulation  of  this  condition  is  given 
in  Section  4.4. 

This  work  is  organized  in  4  Sections  and  3  Appendices.  In  Section  §1.  we  recall  some 
basic  results  about  Hamiltonian  matrices.  In  particular  we  discuss  the  canonical  form 
for  Hamiltonian  matrices.  In  Section  §2  we  define  the  concept  of  eigenvalue  of  positive, 
negative  and  indefinite  type,  and  we  study  some  consequences  related  to  a  perturbed 
eigenvalue  problem  for  the  canonical  blocks.  In  Section  §3  we  study  a  generalized 
eigenvalue  problem  associated  to  (0.3)  and  we  prove  a  proposition  fundamental  in  the 
proof  of  Theorem  I.  Then  we  prove  Theorem  I.  In  section  §4  we  describe  the  reduction 
to  the  center  manifold  and  we  prove  Theorem  II. 

In  Appendix  A  we  describe  some  basic  properties  of  the  logarithm  of  matrices.  In 
Appendix  B  we  go  briefly  through  perturbation  theory  for  eigenvalue  problems  with 
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symmetric  matrices  and  we  apply  it  to  give  a  proof  to  the  results  in  Section  §3.  And 
in  Appendix  C  we  present  the  changes  necessary  to  the  proof  of  the  Center  Manifold 
Theorem  in  order  to  suit  our  needs. 


I  want  to  express  my  gratitude  to  my  thesis  advisor  Professor  Paul  Rabinowitz  for 
his  guidance,  encouragement  and  patience.  Also  I  want  to  thank  professor  E.  Zehn- 
der  who  calls  to  my  attention  the  paper  by  Laub  and  Meyer  on  canonical  form  for 
Hamiltoniem  matrices. 
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§l.BACKGROUND  ON  HAMILTONIAN  MATRICES 

In  this  section  we  will  describe  some  basic  facts  about  Hamiltonian  matrices.  In 
particular  we  will  present  the  canonical  form  of  a  Hamiltonian  matrix;  this  canonical 
form  show  the  symplectic  invariants  we  will  use  in  our  study  of  subharmonics  solutions 
for  (HS). 

A  constant  real  matrix  C  satisfying 

JC^J=-C  (1.1) 

is  called  Hamiltonian  matrix.  Here  J,  as  usual,  denotes  the  standard  symplectic  struc¬ 
ture  in 

■'=(-/  O' 

and  r  denotes  the  transpose  of  a  matrix.  We  note  that  if  B  is  a  symmetric  matrix,  then 
JjB  is  a  Hamiltonian  matrix,  and  conversely,  if  C  is  a  Hamiltonian  matrix  then  B  = 
JC  is  symmetric.  Hamiltonian  matrices  correspond  to  the  linear  part  of  autonomous 
Hamiltonian  systems. 

If  C  is  the  matrix  of  a  linear  system 

i  =  Cz  (1.3) 

where  C  is  Hamiltonian,  then  we  can  transform  (1.3)  through  a  linear  change  of  coor¬ 
dinates  Q  without  changing  the  Hamiltonian  structure  if  Q  satisfies 

Q^JQ  =  J.  (1.4) 

A  matrix  Q  satisfying  (1.4)  is  called  a  symplectic  matrix. 

The  transformation  Q~^CQ,  with  Q  a  symplectic  matrix,  defines  a  simileirity  rela¬ 
tion.  In  such  a  situation  canonical  forms  are  of  great  interest. 

Williamson,  in  [15],  obtained  the  bcisic  results  for  canoniced  forms  of  Hamiltonian 
matrices  .  A  more  constructive  presentation  of  canonical  forms  is  proposed  by  Laub 
and  Meyer  [10].  We  will  tahe  the  canonical  form  as  presented  in  the  paper  of  Bryuno 
[4],  there  simpler  canonical  forms  than  in  [10]  are  given. 

Before  discussing  the  canonical  forms,  we  state  some  basic  lemmas  on  Hamiltonian 
matrices.  For  the  proof  we  refer  the  reader  to  [10], 

Lemma  1.1. 

If  X  is  an  eigenvalue  of  the  Hamiltonian  matrix  C  then  —A,  A  and  —A  are  also 
eigenvalues  of  C  with  the  same  multiplicity  as  A.  In  particular  0  has  even  multiplicity 
if  it  is  an  eigenvalue  of  C. 

With  respect  to  the  eigenspace  associated  to  an  eigenvalue  of  a  Hamiltonian  matrix 
we  have  the  following  lemma. 
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Let  us  define 


(1.5) 

(1.6) 


ry*.(A)  =  ker(C  —  AJ)*,  and 

nW  =  U 

k>i 

Lemma  1.2. 

If  C  is  a  Hamiltonian  matrix  and  A  is  an  eigenvalue  of  C  then  if  fi  is  another 
eigenvalue  of  C,  fi  ^  —A  implies 

(z,  Jw)  =  0  Vz  6  J7(A),  'iw  E  77(A). 


Corollary  1.1. 

If  C  is  a  Hamiltonian  matrix  and  Q  is  an  invertible  matrix  so  that 

J  =  Q-'CQ 


(1.7) 


is  the  real  Jordan  canonical  form  of  C,  then 

Q^JQ  =  diag{qi,q2,...,qp)  (1.8) 

is  a  block  diagonal  matrix  with  one  block  qi  for  each  pair  (A,  —A)  of  eigenvalues  of  C . 

Even  though  the  real  Jordan  canonical  form  of  a  Hamiltonian  matrix  already  shows 
some  special  structure,  it  is  not  enough  for  our  purposes.  The  next  theorem  gives 
the  (real)  canonical  form  of  a  Hamiltonian  matrix  when  only  symplectic  changes  of 
coordinates  are  allowed.  The  proof  can  be  found  in  [4]  and  [15]. 

The  following  notation  will  be  useful  in  this  and  the  next  sections. 


/o  .. 

,.  0 

1  1 

/o  ., 

..  0 

1\ 

1  .. 

,.  0 

0 

0  ., 

,.  1 

0 

; 

; 

II 

: 

; 

•  0 

,.  1 

0) 

1  ' 

,.  0 

0/ 

d  = 


A  =  diag(l,0,  ...,0)  and  I  =  !„■ 

The  size  of  each  matrix  will  be  easily  determined  by  the  context. 


(1.9) 


(1.10) 


Theorem  1.1. 

If  C  is  a  Hamiltonian  matrix  then  there  exists  a  symplectic  matrix  Q  such  that 


Q-^CQ  = 


/Dll 

0 

...  0 

D21 

0 

...  0  \ 

0 

Di2 

...  0 

0 

D22 

...  0 

0 

0 

...  Du 

0 

0 

. . .  D23 

Dzi 

0 

...  0 

Dai 

0 

0 

0 

D32 

...  0 

0 

Da2 

0 

\  0 

0 

...  Du 

0 

0 

...  Das/ 

(1.11) 
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where  the  matrix  formed  by  combining  blocks  in  the  same  relative  position 

jy  fDu  £>2  A 

‘  V^3.  £>4,; 


(1.12) 


has  a  simple  form  analogous  to  the  real  Jordan  canonical  blocks.  Every  block  is  as¬ 
sociated  to  one  eigenvalue  of  C,  but  there  can  be  more  than,  one  block  for  a  speciBc 
eigenvalue.  More  precisely  the  matrices  Di,  that  we  call  generically  D,  can  have  one  of 
the  following  forms: 

1.  Block  of  Cl  type. 

Block  associated  to  an  eigenvalue  A  =  a  + 1/?,  ft  >  0,  /?  >  0  of  multiplicity  1. 


D  = 


^al  +  d  —^I 
pi  al-l-d 
0  0 

\  0  0 


(1.13) 


Here  D  is  a  4l  x  41  matrix. 

2.  Block  of  C2  type. 

Block  associated  to  an  eigenvalue  X  =  a,a  >  0  of  multiplicity  1. 


D 


al  d  0 

0  — ftJ  —  d^ 


(1.14) 


Here  D  is  a  21  x  21  matrix. 

3.  Block  of  C3  type. 

Block  associated  to  an  eigenvalue  X  =  i/?,  P  >  0  of  multiplicity  1. 


D  =  a 


0  K{pl  +  d^) 

-K{pl  +  rf)  0 


(1.15) 


Here  D  is  a  21  x  21  matrix.  The  number  a  can  be  1  or  —1,  and  it  is  an  additional 
symplectic  invariant. 

4.  Block  of  C4  type. 

Block  associated  to  the  eigenvalue  A  =  0  of  multiplicity  1  odd. 


D  = 


d  0 
0  -d'- 


(1.16) 


Here  D  is  a  21  x  21  matrix. 

5.  Block  of  C5  type. 

Block  associated  to  the  eigenvalue  A  =  0  of  multiplicity  1. 


£>  = 


d  a  A 
0  -d^ 


(1.17) 
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Here  D  is  a  21  x  21  matrix.  The  number  a  can  be  1  or  —1,  and  it  is  an  additional 
invariant. 

Remark  1.1.  If  we  compare  with  the  real  Jordan  canonical  blocks,  we  see  that  the 
canonical  blocks  described  in  the  theorem  are  more  complicated  and  cases  C3  and  C5 
carry  an  extra  invariant,  namely  a. 

In  what  follows  we  only  be  interested  in  blocks  of  C3,  C4,  and  C5  type,  i.e.  in 
blocks  corresponding  to  purely  imaginary  eigenvalues. 

The  following  lemmas  are  rather  obvious,  but  we  state  them  for  later  reference. 

In  the  case  of  a  block  of  C3  type  we  can  write 


where 


D  =  -a^V  -  aW 


and 


(1.18) 

(1.19) 


Lemma  1.3. 

1.  V  and  W  are  Hamiltonian  matrices. 

2.  V  J  and  W  J  are  symmetric  matrices. 

3.  V  commutes  with  W  and  J. 

4.  V  is  invertible,  and  W  is  nilpotent. 

Proof.  Parts  1.,  2.  and  4.  are  direct  and  part  3.  requires  an  easy  computation. 
Lemma  1.4. 


exp(27rV')  =  I  (1-20) 

Proof.  Noting  that  KK  =  I,  for  a  G  R  we  have 

,  ,,,  ,  aV  a^I 

exp{aV)  =  - 2r 

=  cos(a)/  +  sin(Q)K<I> 


§2.  PERTURBED  EIGENVALUE  PROBLEM  FOR 
THE  CANONICAL  BLOCKS. 

This  section  is  devoted  to  the  definition  of  eigenvalue  of  positive,  negative  or  indef¬ 
inite  type,  in  the  case  of  purely  imaginary  eigenvalues  of  Hamiltonian  matrices. 

The  definitions  we  are  going  to  make  are  motivated  by  the  application  to  subhar¬ 
monic  solutions  of  Hamiltonian  systems  that  we  will  present  in  next  section.  However, 
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we  can  already  see  the  meaning  of  these  definitions  in  studying  some  perturbed  eigen¬ 
value  problems  for  the  canonical  blocks  .  In  Section  §3  we  will  see  how  this  is  related 
to  the  existence  of  subharmonics  for  the  Hamiltonian  system  (HS). 

§2.1.  Definitions. 

In  the  following  definitions  we  assume  we  have  a  Hamiltonian  matrix  C,  and  we 
have  obtained  a  symplectic  transformation  Q  that  taJces  C  into  its  canonical  form.  D 
will  denote  a  generic  block  in  the  canonical  form  of  C  as  given  by  Theorem  1.1. 

Definition  2.1.  We  say  that  a2lx2l  block  D  is  positive  (negative)  if  one  of  the  following 
is  true: 

1.  D  is  a  block  of  C3  type,  1  is  even  and  <7  =  1  (cr  =  —1). 

2.  D  is  a  block  of  C5  type  and  a  =  1  (cr  =  —1). 

Definition  2.2.1/  A  is  a  purely  imaginary  eigenvalue  of  the  Hamiltonian  matrix  C  we 
say  that  X  is  of  positive  type  (negative  type)  if  all  blocks  associated  to  X  are  positive 
(negative).  We  say  that  X  is  of  indehnite  type  if  X  is  not  of  positive  type  nor  of  negative 
type. 

Remark  2.1.  The  concepts  of  positive,  negative  and  indefinite  type  are  symplectic 
invariants. 

Remark  2.2.  See  Remark  3.6  for  a  perturbation  result  for  blocks  of  positive  type  and 
of  negative  type. 


§2.2.  Perturbed  eigenvalue  problems. 

Given  a  block  D  of  the  canonical  form  of  a  Hamiltonian  matrix,  that  is  of  C3,  C4 
or  C5  type  we  will  define  an  associated  eigenvalue  problem  depending  on  a  parameter  e. 
Using  perturbation  theory  we  can  study  the  behaviour  of  the  eigenvalues  of  the  problem 
as  a  function  of  e.  In  particular  we  are  interested  in  the  sign  of  the  small  eigenvalues 
for  small  values  of  the  parameter  e.  The  results  are  given  in  Propositions  2.1,  2.2  and 
2.3.  We  delay  the  proof  of  these  propositions  to  Appendix  B. 

Case  1.  Block  of  C3  type.  Using  the  notation  given  in  Section  §1  we  can  write  a 
block  of  C3  type  in  the  form 

D  =  -adV  -aW  (2.1) 

where  V  and  W  were  defined  in  (1.19). 

Let  us  consider  the  following  eigenvalue  problem 

—[cV  -t-  (tW)v  -f-  XJv  =  0  (2.2) 

where  c  €  R  and  v  €  We  study  the  eigenvalue  A  as  a  function  of  e. 
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This  problem  is  equivalent  to 


{eV  J aW  J)w  +  \w  =  0,  w  =  Jv.  (2.3) 

From  Lemma  1.3  we  know  that  this  is  a  symmetric  eigenvalue  problem,  so  that  the 
behaviour  of  the  eigenvalues  of  (2.3)  nezir  e  =  0  can  be  determined  by  using  the  pertur¬ 
bation  theory  for  such  problems. 

The  simplest  case  occurs  when  the  block  D  in  question  is  2  x  2  ,  i.e.  1=1.  Then 
we  have 

VJ=(^l  1)  and  “). 

Obviously  0  is  cin  eigenvalue  oi  W  J  of  multiplicity  2,  and  for  e  ^  0,  (2.3)  has 

A(e)  =  -£  (2.4) 


as  a  double  eigenvalue. 

The  situation  in  the  general  case  is  more  complicated.  It  is  easy  to  see  that  the 
matrix  W  J  has  0  as  an  eigenvalue  of  multiplicity  2.  We  want  to  determine  what  happens 
to  this  eigenvalue  for  small  e  ^  Q.  From  Theorem  B.l  we  know  that  the  eigenvalues  of 
(2.3)  are  analytic  functions  of  c  near  e  =  0.  Thus  A  =  0  will  be  continued  as 

Aji(e)  =  Aje  -f-  A^e^  4- ...  j  =  1,2.  (2.5) 

Since  the  matrix  V  is  invertible,  W  is  nilpotent  and,  V  and  W  commute,  the  matrix 
(eV  4-  aW)J  is  invertible.  Consequently  in  the  series  (2.5)  there  is  at  least  one  nonzero 
coefficient.  The  following  proposition  determine  the  first  nonzero  coefficient. 

Proposition  2.1. 

If  the  block  D  is  21  x  21  then  the  series  (2.5)  can  be  written  as 

kj(e)  =  <T'-^(-ey  +  ...  j  =  l,2.  (2.6) 


Proof.  See  Appendix  B.O 

Remark  2.3.  If  is  a  positive  block  of  C3  type  then  according  to  Definition  2.1  /  is 
even  and  a  =  1,  this  implies  that 

Aj(€)  =  e' 4- ...  i  =  l,2 

is  positive  for  e  7^  0  in  a  neighborhood  of  0.  If  D  is  a  negative  block  then  Ay(e)  is 
negative  near  zero.  If  D  is  not  positive  nor  of  negative  type  then  Aj(e)  takes  positive 
and  negative  values  near  0. 


Case  2.  Block  of  C4  type.  We  recall  that  a  block  of  C4  type,  according  to  Theorem 
1.1  has  the  form 


here  D  is  a  2/  x  21  matrix,  with  I  being  odd. 

As  in  the  case  just  treated  we  will  consider  an  eigenvalue  problem 

(eil  +  D)v  +  XJv  =  0,  (2.8) 

with  V  €  and  e  G  R,  or  equivalently 

(—eiJ  —  DJ)w  +  Xw  =  0,  w  =  Jv.  (2-9) 

The  matrix  DJ  is  easily  seen  to  be  symmetric,  and  iJ  is  hermitian.  We  can  treat  (2.9) 
using  perturbation  theory  as  before.  The  eigenvalues  continuing  A  =  0  have  the  form 
(2.5)  and  the  next  proposition  determine  the  first  nonzero  coefficient  in  the  series. 

Proposition  2.2. 

If  the  block  D  is  21  x  21  of  C4  type  then  the  series  (2.5)  can  be  written  as 

Xj(e)  =  (-iye‘  +  ...  j  =  l,2.  (2.10) 


Proof.  See  Appendix  B.  <C> 

Remark  2.4.  We  observe  that  for  e  small  Aj(e),  j  =  1,2  takes  positive  and  negative 
values.  According  to  Definition  2.1  a  block  of  C4  type  is  not  positive  nor  negative. 

Case  3.  Block  of  C5  type.  Now  we  turn  to  study  blocks  of  the  C5  type,  where  again 
the  invariant  cr  appears. 

From  Theorem  1.1,  a  block  of  this  type  is  of  the  form 


D  = 


d  crA  'N 
0  -d^ ) 


with  D  a  21  X  2/  matrix  and  cr  =  1  or  —1. 

We  consider  the  eigenvalue  problem 

(ei/ +  £>)?;  + A  Jv  =  0,  (2.11) 

with  V  e  R^^  and  c  G  R,  or  equivalently 

{—tiJ  —  DJ)w  +  An;  =  0,  w  =  Jv.  (2.12) 

DJ  is  symmetric  and  iJ  is  hermitian.  We  apply  again  perturbation  theory. 
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Proposition  2.3 

If  the  block  D  is  21  x  21  of  C5  type,  the  eigenvalue  of  (2.12)  continuing  A  =  0,  can 
be  written  as 

X(e)  =  ^  ^2.13) 

Proof.  The  proof  of  this  proposition  goes  through  the  same  lines  as  as  the  proof  of 
Proposition  2.1  and  2.2.  We  omit  it.O 

Remark  2.5.  By  analyzing  (2.13)  in  the  case  Z)  is  a  positive  block  we  can  see  that  A(e) 
is  positive  in  a  neighborhood  of  e  =  0.  Similarly  if  D  is  negative  then  A(e)  is  negative 
near  c  =  0. 

§2.3  ivr  as  a  special  eigenvalue  of  a  Hamiltonian  matrix 
If  we  want  to  study  a  Hamiltonian  system  like 

i  =  JA(t)z  +o(|2l), 

and  we  assume  that  A(t)  =  X  is  constant  then  J A  is  a  Hamiltonian  matrix  and  we  can 
apply  the  Definitions  2.1  and  2.2  to  JA.  If  A{t)  is  time  dependent  then,  using  Floquet 
theory  we  can  transform  the  system  to  one  with  a  constant  linear  term.  However  if  the 
system  has  negative  Floquet  multipliers  we  need  to  deal  with  a  new  situation. 

Hamiltonian  matrices  can  be  obtained  from  symplectic  matrices  by  taking  the  log¬ 
arithm.  A  special  situation  arises  when  the  symplectic  matrix  has  negative  eigenvalues. 

To  begin  we  recall  some  results  that  are  the  symplectic  analogues  of  Lemmas  1.1 
and  1.2.  The  proofs  can  be  found  in  [10]. 

Lemma  2.1. 

If  A  is  an  eigenvalue  of  a  symplectic  matrix  X  then  1/A,  A,  1/A  are  also  eigenvalues 
of  X  with  the  same  multiplicity  as  X.  In  particular  1  and  —1  have  even  multiplicity  if 
they  are  eigenvalues  of  X . 

Let  us  define 

ik{X)  =  ker(X  -  A/)^  and  (2.14) 

«A)  =  U  6(A),  (2.15) 

ik>l 

and  similarly 

^jt(A)  =  ker(J\r“*  —  A/)^,  and  (2.16) 

m  =  U  (2-17) 

*>1 

Then  we  have 
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Lemma  2.2. 

If  X  is  a  symplectic  matrix  and  X  is  an  eigenvalue  of  X  then 

1.  e(A)  =  |(l/A) 

2.  If  y.  is  another  eigenvalue  of  X  then  y  ^l/\  implies 

(z,Jw)  =  0  Vz  e  ^(A),  '^w  e  ({y). 


Corollary  2.1. 

If  X  is  a  symplectic  matrix  and  Q  is  an  invertible  matrix  so  that 

J  =  Q~^XQ  (2.18) 

is  the  real  Jordan  canonical  form  of  X,  then 

=  diag(9i,52,-.,9p)  (2.19) 

is  a  block  diagonal  matrix  with  one  block  qi  for  each  pair  (A,  1/A)  of  eigenvalues  of  X. 
For  X  symplectic  let  us  consider 


C  =  log(X)  (2.20) 

where  log  denotes  the  logarithm  defined  in  Appendix  A.  If  X  does  not  have  eigenvalues 
on  the  negative  real  axis  then  C  is  a  real  Hamiltonian  matrix  (See  Lemma  A. 2  and  A. 3). 
However,  if  X  has  negative  real  eigenvalues  C  is  no  longer  real.  The  matrix  C  is  not 
Hamiltonian  in  the  sense  we  defined  before,  but  it  satisfies 

JC*J  =  C  (2.21) 

where  *  denotes  the  conjugate  transpose  of  C.  For  such  a  matrix  C  we  want  to  have 
‘real’  canonical  forms.  We  will  reduce  this  question  to  Theorem  1.1. 

Let  us  assume  that  X  has  2p  negative  eigenvalues,  and  let  us  reorder  J  as 

^'=(o  ",)  (2.22) 

where  Jp  and  Jg  are  the  standard  symplectic  structures  in  and  respectively. 
Here  q  =  n  —  p.  From  Lemma  2.2  there  is  an  invertible  matrix  Q  so  that 

Q^JQ  =  J  (2.23) 

and 

Y  =  Q-^XQ  (2.24) 
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is  a  block  matrix 


with 


■  y 


(2.25) 


YrjjYj  =  Jj,  j  =  p,q 


(2.26) 


and  Yp  has  only  negative  eigenvalues,  while  Yq  does  not  have  negative  eigenvalues. 
Let 

5j=log(yj),  j  =  p,q 


then,  by  Lemma  A. 3,  we  have 

JjSjJj  —  Sj,  j  =  p,  9. 

Lemma  A.2  guarantees  that  Sq  is  real,  and  Lemma  A.4  gives  that 

Sp  —  tTvJ'^p  "I”  Sp 

where  Sp  =  log(— 1^)  is  a  real  matrix,  and  l2p  is  the  identity  in 
From  (2.28)  we  see  that  Sp  satisfies 


JpSpJp  —  Sp- 


If  we  define 


and  5  = 


then  we  have 

and  the  real  matrix  S  satisfies 

If  we  define 


■  ('? :) 

log(y')  =  i-Ki  +  s 

JS^J  =  s. 


Sp  0 
0 


(2.27) 


(2.28) 


(2.29) 


(2.30). 


(2.31) 

(2.32) 


Ci^itQIQ-^  and  Cr  =  QSQ~'^  (2.33) 

then  from  (2.31)  ,  (2.32)  and  the  block  structure  of  /,  J  and  S  the  following  hold 


C  =  iCi  +  Cr,  (2.34) 

JC;j  =  Cr,  (2.35) 

JC^J  =  -Ci,  and  (2.36) 

CiCr  =  CrCi-  (2.37) 
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If  —1  is  an  eigenvalue  of  X  of  multiplicity  2k  then  0  is  an  eigenvalue  of  multiplicity  2k 
of  Sp  and  from  (2.30)  Sp  is  a  Hamiltonian  matrix.  We  have  the  following  extension  of 
Definition  2.2. 


Definition  2.3. Let  X  be  a  symplectic  matrix,  and  C  =  log(Ar).  If  —1  is  an  eigenvalue 
of  X  of  multiplicity  2k,  then  ix  is  an  eigenvalue  of  C  of  multiplicity  2k  and  we  say  that 
ITT  is  of  positive  type  (negative  type  or  indefinite  type)  if  and  only  if  0  is  an  eigenvalue 
of  positive  type  (negative  type  or  indeGnite  type)  of  the  matrix  Sp. 


§3.PROOF  OF  THEOREM  I. 

In  this  section  we  begin  the  study  of  the  Hamiltonian  system  (HS) 

i  =  JA(t)z  +  JHz{z,t) 

where  the  Hamiltonian  function  is 

H(z,t)  =  ^iA(t)z,z)  +  Hiz,t)  (3.1) 

and  we  assume  it  satisfies  the  basic  hypotheses; 

(HO)  H  is  of  class  near  z  =  0, 

(HI)  H(0,  #)  =  0  and  /f^(0,  t)  =  0  V<  €  R, 

(H2)  H  is  1-periodic  in  the  t-variable. 

Assuming  that  H  is  1-periodic  instead  of  T-periodic  does  not  restrict  the  analysis. 
Under  these  hypotheses  (HS)  represents  the  Hamiltonian  system 

i  =  JHz{z,t),  (3.2) 

A(t)  is  the  Hessian  matrix  of  H  with  respect  to  z  at  z  =  0,  and  H  represents  the  higher 
order  terms  of  H,  H(z,t)  =  o(\z\^)  uniformly  in  t,  for  z  near  0. 

Let  us  consider  the  linear  system 

z  =  J  A{t)z.  (3.3) 

If  X{i)  is  the  fundamental  matrix  if  system  (3.3),  i.e. 

X{0)  =  hn  (3.4) 

X{t)  =  JA{t)X{t)  (3.5) 

we  define  X  =  X{1)  and  C  =  log(X),  where  the  logarithm  is  defined  in  Appendix  A. 
The  eigenvalues  of  C  are  called  the  Floquet  exponents  of  (3.3).  We  will  assume 
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(H3)  There  is  at  least  one  Floquet  exponent  of  (3.3)  on  the  imaginary  axis  that  is 
not  of  negative  type. 

In  other  words,  (H3)  says  that  C  has  at  least  one  purely  imaginary  eigenvalue  that  is 
not  of  negative  type.  Finally  we  assume 

(H4)  There  are  constants  c>  Q,r  >  +  2  >  p  >  q  >  2  such  that 

{H,{z,t\z)>qH{2,t)>c\zY  (3.6) 

Vr  e  \z\  <  r,  Vt  e  R. 

In  this  section  we  will  prove 
Theorem  I 

If  (HO),  (HI),  (H2),  (H3)  and  (H4)  hold,  then  there  is  an  integer  N  so  that  for 
every  k  >  N,  there  exist  a  nontrivial  solution  Zf,  of  (HS)  with  period  k,  and  {zik}it>iv 
converges  to  zero  as  fc  — >  oo  uniformly  in  C^(R,R^”). 

Remark  3.1.  In  Section  §4  we  will  see  that  adding  some  differentiability  to  H  we  will 
be  able  to  prove  the  result  of  Theorem  I  under  less  restrictive  conditions  than  (H4), 
namely  we  need  only  to  assume  (3.6)  holds  on  the  center  manifold  of  (HS). 

The  proof  of  Theorem  I  is  presented  in  the  next  two  sections.  We  begin  in  Section 
§3.1  by  proving  a  result  on  a  generalized  eigenvalue  problem  associated  to  (3.3).  For 
that  purpose  we  transform  (3.3)  to  a  linear  system  with  constant  coefficients  and  then 
we  analize  the  resulting  system  in  terms  of  the  canonical  form  associated  to  it.  Next  we 
study  the  nonlinear  problem  in  Section  §3.2  using  variational  methods. 


§3.1  The  Linearized  Problem. 

The  proposition  we  will  prove  now  is  the  basis  of  the  proof  of  Theorem  I.  It  is 
here  that  hypothesis  (H3)  is  crucial.  We  will  first  reduce  system  (3.3)  to  a  system 
with  constant  coefficients  via  Floquet  Theory.  Because  the  matrix  X  =  A’(l)  can 
have  negative  eigenvalues  we  will  need  to  introduce  complex  matrices  or  change  our 
periodicity  condition  on  the  solutions.  We  prefer  the  second,  because  in  Section  §4  this 
will  be  more  convenient. 

Let  i?  be  a  matrix  so  that  H?  =  I. 

Definition  3.1. 

Let  k  £N,  a  function  z  :  R  — >  R^"  is  said  to  be  Rk-periodic  if 

z{t  +  k)  =  R*z(t)  Vt  €  R.  (3.7) 


We  note  that  every  R-periodic  function  is  2-periodic. 

Lemma  3.1. 
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There  are  matrices  Cr,i2  and  P(i)  so  that  Cr  is  Hamiltonian,  B?  =  I,  Pit)  sym- 
plectic  and  P(t  +  1)  =  P(<)i2  Vt  6  R  for  which  we  have 
z(t)  is  a  k-periodic  solution  of  (3.3)  if  and  only  if 


is  an  Rk -periodic  solution  of 


^(t)  =  P(t)  z{t) 


i  = 


(3.8) 


Proof.  The  proof  is  an  application  of  Floquet  theory.  Let  Xit)  be  the  fundamental 
matrix  of  (3.3)  as  defined  above.  Since  (3.3)  is  Hamiltonian  ^(t)  satisfies 


X{t)^JX{t)  =  J  V<  6  R. 


(3.9) 


To  obtain  (3.9)  we  differentiate  the  matrix  2(1)  =  X(t)'^JX(t)  and  use  (3.4),  (3.5)  and 
the  fact  that  A(t)  is  symmetric.  Identity  (3.9)  says  that  X(t)  is  symplectic  for  every 
f  €  R,  consequently  X  =  X(l)  is  symplectic. 

By  Lemma  A. 3  there  is  a  matrix  C  so  that 


X  =  exp(C),  and 
C  =  JC*J. 


(3.10) 

(3.11) 


As  we  have  mentioned  before  C  is  not  necessarily  real  because  X  can  have  negative 
eigenvalues.  However  we  can  always  decompose  C  as  C  =  iCi  +  Cr  with  C,-  and  Cr 
satisfying  (2.35),  (2.36)  and  (2.37). 

Let  us  define 

Pit)  =  Xit)expi-Crt) 

then  Pit)  is  invertible  and  obviously 


A'(t)  =  P(t)exp(Crt). 

Since  A’(t)  is  symplectic  and  Cr  satisfies  (2.35)  P(<)  is  also  symplectic,  i.e. 

P(t)’‘/P(t)  =  J. 

On  the  other  hand,  since  (3.3)  is  1-periodic  it  is  easy  to  see  that 

Xit  +  1)  =  Xit)Xil). 

Therefore 

P(<  +  1)  =  Xit  +  l)exp(-Cr(t  +  1)) 

=  Xit)X  exp(— C)exp(— Cr<  +  iCi) 

=  A'(t)exp(— (7rt)exp(iCi) 

=  P(<)exp(fCi). 


(3.12) 

(3.13) 


(3.14) 


(3.15) 
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where  we  used  (2.37)  and  the  definition  of  P{t)-  If  we  define 

iZ  =  exp(iC'i)  (3.16) 

then  from  the  definition  of  C,-,  (2.33), 

=  /  (3.17) 

and  what  we  showed  in  (3.15)  becomes 

P{t+1)  =  P{t)R.  (3.18) 

Finally  we  have  the  following  identity  obtained  from  differentiating  (3.14)  and  using 

(3.5) 

P(t)  +  P(t)Cr  =  JA{t)P{t).  (3.19) 

Let  z{t)  be  a  /'-periodic  solution  of  (3.3),  and  ^(t)  =  P{t)~^  z{t).  Then  by  differentiating 
z(t)  =  P(t)^(t)  and  using  (3.3)  we  obtain 

p(l)«t)  +  P(t)m  =  J>l(t)P(i)«i).  (3.20) 

Using  (3.19)  in  (3.20)  and  multiplying  by  P(t)~^  we  have 

e  =  Cr^- 

Since  P{t)  satisfies  (3.18),  and  z  is  ^-periodic  we  have 

«t+i)  =  ij*e((). 

The  reverse  implication  is  also  easily  obtained.^ 

Next  we  prove  a  lemma  that  will  be  used  in  Section  §4,  which  is  a  sort  of  inverse 
of  Lemma  3.1  in  the  sense  that  finding  /Zfc-periodic  solutions  of  a  system  like  (3.8)  is 
equivalent  to  find  /;-periodic  solutions  of  a  system  with  variable  coefficients.  Note  that 
this  is  not  given  by  Lemma  3.1  because  in  this  situation  we  do  not  have  the  fundamental 
matrix  X{t)  available. 

We  defined  iZ  in  (3.16)  in  terms  of  C,.  We  recall  that  C,-  was  defined  in  (2.33)  eis 

Ci  = 

where  Q  is  a  symplectic  matrix,  and  /  is  given  by 


and  l2p  is  the  identity  in 
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We  can  write  the  matrix  i?  in  a  different  way,  and  we  do  it  now.  If 


11 

(3.21) 

with 

(3.22) 

then  it  is  easy  to  see  that  for  the  matrix 

Ci  =  irQIQ-^ 

(3.23) 

we  have 

R  =  exp(C'i). 

(3.24) 

Lemma  3.2. 

^(t)  is  an  Rk-periodic  solution  of 


if  and  only  if  w(t)  =  R(t)^(t)  is  a  k-periodic  solution  of 

w  =  (R(t)  +  i?(<)Cr)i2~^(t)u; 


where  R(t)  =  exp(Cif). 

Proof.  By  direct  calculations.^ 

Remark  3.2.  We  note  that  in  general  the  matrices  Ci  and  Cr  do  not  commute.  That 
is  the  reason  why  we  cannot  use  Cr  +  Ci  as  the  logarithm  of  X. 

The  following  proposition  is  fundamental. 

Proposition  3.1. 

Assume  the  linear  system  (3.3)  satisfies  (H3).  Then  there  is  a  1-periodic,  continuous 
matrix  function  U(i),  that  is  symmetric  and  positive  dehnite,  such  that  for  every  k  gN 
there  is  a  k-periodic  function  Zk,  and  a  number  Xk,  A*  >  0  and  Xk  =  0{\/k)  as  fc  — >  oo, 
satisfying 

z  =  JA{t)z  +  XJU{t)z.  (3.25) 

In  other  words  the  generalized  eigenvalue  problem  (3.25)  has  at  least  one  eigenvalue 
Ajt  >  0  with  Xk  =  0(1/ k)  as  k  —*  oo. 

Proof.  From  Lemma  3.1  we  can  transform  (3.3)  into 
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with  Cr  a  Hamiltonian  matrix.  Let  N  be  the  canonical  form  of  Cr  as  in  Theorem  1.1. 
Then  there  is  a  Q  so  that 

Q'^JQ  =  J  and  (3.26) 

Q-^CrQ  =  N.  (3.27) 

From  hypothesis  (H3)  we  know  that  N  has  at  least  one  eigenvalue  A  that  is  purely 
imaginary  and  not  of  negative  type.  For  that  eigenvalue  there  exists  a.  21  x  21  block  D 
that  is  not  negative. 

Let  /b  €  N,  and  consider  the  following  eigenvalue  problem 

X  =  Dx  +  A/r,  x(0)  =  x{k)  (3.28) 

where  x  £  r2',a  g  r  and  J  is  the  standard  symplectic  structure  in  R^L  We  do  not 
distinguish  the  dimension;  it  will  be  clear  from  the  context. 

The  block  D  can  be  of  C3,  C4,  or  C5  type.  We  analyze  each  case  separately. 
CASE  1.  Block  of  C3  type. 

In  this  case  A  =  i/3,/3  >  0,  and  the  matrix  D  has  the  form  given  by  (1.18) 


D  =  -al3V  -  aW  (3.29) 

Let  a  €  [— 27r,  27r]  so  that  exp(—ia)  =  exp(—icrjdk);  we  determine  the  sign  of  a  later. 
The  eigenvalue  problem  (3.28)  is  related  to 

— (qV  +  kaW)v  +  XkJv  =  0.  (3.30) 

In  fact,  let  us  assume  (A,u)  satisfies  (3.30)  and  let  us  define 

Xk{t)  =  exp((— —  aW  +  \J)i)  v  (3.31) 

Certainly  ijt(t)  satisfies  the  differential  equation  in  (3.28),  but  it  also  satisfies  the  bound¬ 
ary  condition 

Xk{k)  =  exp(fc(— cr^V  —  crW  -|-  AJ))  v 

=  exp(—crj3kV)exp(—akW  -1-  XkJ)v 
=  exp(— aF)exp(— crArW  -|-  XkJ)v 
=  exp(— qV  —  akW  XkJ)v 
=  u, 

here  we  used  Lemmas  1.3  and  1.4,  the  definition  of  a  and  (3.30).  Thus  it  is  enough  to 
study  (3.30). 

Dividing  (3.30)  by  k  we  get 


-i^V  +  aW)v  +  XJv  =  0 

K 


(3.32) 
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that  is  exactly  the  eigenvalue  problem  (2.2),  with  e  =  ajk.  Because  D  is  not  negative, 
considering  a  with  the  adequate  sign,  from  Proposition  2.1  we  have  that  for  every  k 
laxge  enough  there  is  a  Ajt  >  0,  an  eigenvalue  of  (3.32)  and  hence  of  (3.30).  Moreover 
Ajt  =  0{l/k). 

Thus,  for  every  G  N  large  enough,  (3.28)  has  an  eigenvalue  At  >  0,  with  Ajt  = 
0(l/k)  and  associated  eigenfunction  Xk  that  is  A:-periodic. 

The  remaining  cases  correspond  to  blocks  of  C4  and  C5  type.  These  blocks  have 
two  possible  origins,  they  can  come  from  an  eigenvalue  0  of  C  or  from  the  eigenvalue  ztt 
of  C.  We  distinguish  between  the  two  possibilities. 

CASE  2.  Block  of  C4  type  (  A  =  0  is  an  eigenvalue  of  C)  In  this  case  the  matrix  D 
has  the  form 


Let  a  =  27r  and  consider  the  eigenvalue  problem 

(ail  +  kD)v  -f  \kJv  =  0.  (3.33) 

If  (A,t;)  is  a  solution  of  (3.33)  then  the  function 

Xk(t)  =  exp((Z?  +  XJ)t)  V  (3.34) 

satisfies  equation  (3.28).  Using  (3.33)  we  see  that 

Xk(k)  =  exp{kD  -f  XkJ)v 

=  exp(27ri/  +  kD  +  XkJ)  v 

=  V 

so  that  Xk  is  fc-periodic.  Thus  it  is  enough  to  solve  (3.33),  but  dividing  (3.33)  by  k  we 
see  that  it  corresponds  exactly  to  the  eigenvalue  problem  (2.8)  with  e  =  2tz fk 

Stt 

{^iI  +  D)v  +  XJv  =  0 

hf 

and  Proposition  2.2  guarantees  that  for  large  k  there  is  a  At  >  0,  Afc  =  0(1/ k),  an 
eigenvalue  of  (3.33).  Thus  we  have  obtained  the  same  conclusion  as  in  Case  1. 

CASE  3.  Block  of  C4  type  (corresponding  to  the  eigenvalue  A  =  iir  of  C).  Here  we 
proceed  as  in  Case  2  but  we  will  take 

{TT  if  is  odd 
27r  if  k  is  even 

The  only  difference  with  Case  2  is  that  the  corresponding  function  Xk  will  satisfies 


^k(k)  =  (-l)*a:t(0). 
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CASES  4  and  5.  D  is  a  block  of  C5  type  (corresponding  to  A  =  0  or  itt  of  C). 
This  cases  are  similar  to  Cases  2  and  3  respectively.  We  note  that  D  is  of  positive  type 
so  that  from  Proposition  2.3  we  can  guarantee  that  the  eigenvalue  Xk  can  be  chosen 
positive  and  also  of  order  0(l/fc). 

Now  we  summarize  what  we  have.  From  (2.33)  and  (3.16)  we  see  that 

R  =  QRiQ-^ 

where 


Then,  what  we  have  shown  is  that  the  eigenvalue  problem 

i  =  N^  +  XJ(,  m  =  RH(0)  (3.35) 

has  an  eigenvalue  A/t  >  0,  Ajt  =  0{l/k)  with  eigenfunction  (k,  for  any  ^  G  N  large 
enough. 

Changing  coordinates,  Wk(t)  =  Q^k(i)  satisfies 

w  =  CrW  +  XQJQ~^w,  w(k)  =  R^w{0).  (3.36) 

Finally  defining 

Zkit)  =  P{t)wk{t)  (3.37) 

and 

U{t)  =  ~JP{t)QJQ-^P{t)-\  ■  (3.38) 

using  the  same  argument  of  the  proof  of  Lemma  3.1,  we  see  that 

Zk  =  JA{t)zk  +  XkJU{t)zk,  and 

Zk{k)  =  rjt(O). 

It  is  only  left  to  prove  that  U{t)  is  symmetric  and  positive  definite.  Because  Q  and  P{t) 
are  symplectic  we  have 


Q^J  =  JQ-^  and  P{ty  J  =  J  P{t)-^  (3.39) 

and  then  from  (3.38)  and  (3.39)  we  obtain 

U(t)  =  rP(t)QQ^P(tyj,  (3.40) 

from  where  U{t)  is  clearly  symmetric  and  positive  definite.  A  simple  computation  using 
(2.33),  (3.16)  and  (3.18)  will  show  that  U{t)  is  also  1-periodic.O 
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§3.2  The  Nonlinear  Problem. 

Now  we  give  a  proof  of  Theorem  1.  We  start  with  some  preliminaries. 

Let  us  consider  the  eigenvalue  problem 

—  Jz  —  A(t)z  =  \U{t)z  (3-41) 

in  L^([0,  fc],  R^")  =  L|  with  periodic  boundary  conditions.  Here  U  is  the  matrix  function 
defined  in  Proposition  3.1.  In  Lj.  we  will  consider  two  inner  products 


k 


(u,u)2  =  /  {u(t)^v(t))dt,  and 

Jo 

(3.42) 

<u,v>2=  f  (U{t)u{t),v{t))dt. 

Jo 

(3.43) 

Let  I  •  I2  and  ||  •  II2  be  the  norms  induced  by  (•,-)2  and  <  •,•  >2  respectively. 

Since  the  matrix  function  U(t)  is  periodic,  continuous  and  positive  definite  there  is 
a  constant  a  >  0  such  that 

“  ^1  U(t)  \<  a  Vt  6  R,  and  (3-44) 

a 

-<\U{t)-^\<a  VteR.  (3.45) 

a 

where  |  •  |  denotes  the  L(R^",R^")  norm  induced  on  matrices  (using  the  usual  norm  in 
R^"). 

Lemma  3.3. 

There  are  constants  ci  >  0  and  C2  >  0  independent  of  k  so  that 

I  u  |2<  Cl  11  u  II2  and 
il  u  Il2<  C2  1  u  I2  Vu  G  Ll- 
Proof.  ^ 

I  w  I2  =  /  {u{t),u{t))dt 

Jo 

{U{t)Uit)-'^uit),u{t))dt  (3.48) 

=  <  U{t)~^u,u  >2 

<11  C/-^ull2|!ull2. 


(3.46) 

(3.47) 
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But  from  (3.45) 


\\U  II2  =  /  iu{t),U(t)  ^u{t))dt 

JQ 

<  a  I  u  I2  . 

Thus,  from  (3.48)  and  (3.49)  we  obtain 

\u  \2<  \/a\\u  II2  . 


(3.49) 


On  the  other  hand  ^ 

II  ’x  111  =  /  (U(t)u{i),u(t))dt 

Ja 

<1  Uu  I2I  u  I2  . 

But  ^ 

lUu\l=  [  {U(t)u(t),V{i)u(t))dt 
Jo 

<a^\u  II 

so,  from  (3.50)  and  (3.51)  we  obtain 

11  ti  Il2<  \/a  I  u  I2  . 


(3.50) 


(3.51) 


Taking  ci  =  C2  =  \/a  we  get  (3.46)  and  (3.47). 

For  the  eigenvalue  problem  (3.41)  we  have  the  following  standard  result. 

Lemma  3.4. 

The  eigenvalue  problem  (3.41)  possesses  a  sequence  of  eigenvalues  that 

extends  from— 00  to-foo.  Each  eigenvalue  is  isolated  and  has  £nite  multiplicity.  As¬ 
sociated  to  the  eigenvalues  there  is  a  basis  of  L].,  composed  of  eigenvectors  {vj}j£2, 
orthonormal  with  respect  to  the  inner  product  <,  >2- 

Proof.  We  only  give  the  general  idea  of  the  proof.  First,  by  making  a  transformation 
like  that  of  Proposition  3.1,  and  taking  in  account  the  way  U  was  defined  we  can 
transform  system  (3.41)  into  a  system  with  constant  coefficients.  Second,  by  taking  an 
appropriate  trigonometric  basis  in  we  transform  this  constant  coefficient  differential 
equation  into  a  matrix  one,  where  the  eigenvalues  can  be  calculated.  The  trigonometric 
basis  gives  rise  to  eigenvalues  ±n  with  a  constant  correction  coming  from  the  constant 
term.'^ 

Remark  3.3.  By  Proposition  3.1  there  is  an  eigenvalue  >  0  of  (3.41)  with  A*  = 
0(1/ k)  as  k  — +  00. 

In  order  to  study  the  existence  of  subharmonics  of  (HS)  we  will  study  the  critical 
points  of  a  functional.  Let  E  be  the  Sobolev  space  Wj|!^^’"(R,  R^").  For  details  about 
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the  definition  of  E  and  the  operators  we  define  next  we  refer  the  reader  to  [13].  For 
z  =  (p,  5)  €  and  smooth  we  ;d'^ne 


Q{z)  =  QP-  2 ^(0^)  dt. 


If  C  =  (v’jV’)  ^  E  \s  also  smooth  we  can  define  the  bilinear  form 


B{z,0-  f  ‘<Pp  +  q‘p-{z,A{t)C)dt 

Jo 

=  j  {~Jz~A{t)zX)dt. 

Jo 


Both  Q  and  B  can  be  extended  continuously  to  all  E,  and 


(3.52) 


(3.53) 


^B(z,z)  =  Q{z)  ^zeE.  (3.54) 

Using  Lemma  3.4  we  can  define  the  following  closed  subspaces  of  E 

E'^  =  span{u,7  Vj  >  0} 

E~  =  span{u,/  i/,-  <  0},  and 
E^  =  span{v,7  i/j  =  0}. 

It  is  easy  to  see  that  Q  is  positive  in  E'^,  negative  in  E~  and  zero  in  E°.  Also  for 
z  =  z+  +  2~  +  2*’, 2+  g  E'^,z~  G  E~  and  2®  G  E°  we  have 

Q{z+ +  z- +  z°)  =  Q{z+)  +  Q{z-).  (3.55) 

Moreover,  we  have  the  following  lemma. 

Lemma  3.5. 

If  we  define  ||  2  ||q=  Q{z'^)  —  Q{z~)  +  Jq  \  2°  P  dt  for  every  z  ^  E,  then  there  are 
constants  C3  and  C4  such  that 

C3  11  2  ||q<||  2  ||£;<  C4  II  2  IIq  (3.56) 

i.e.  II  •  IIq  is  an  equivalent  norm  in  E. 

Proof.  For  proving  this  lemma  we  would  proceed  as  in  Lemma  3.4  passing  to  the  equiv¬ 
alent  definitions  of  Q  and  ||  ■  ||q  for  an  equivalent  system  with  constants  coefficients, 
where  the  lemma  is  clear. 0 

As  mentioned  above  in  order  to  find  subharmonics  solutions  of  (HS)  we  will  use  a 
variational  approach.  In  particular  we  will  use  a  version  of  the  generalized  Mountain 
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Pass  Theorem  due  to  Benci  and  Rabinowitz.  We  state  this  theorem  here  and  we  refer 
the  reader  to  [13]  for  its  proof. 

Theorem  3.1. 

Let  E  be  a  real  Hilbert  space  with  E  =  Ei  ®  E2  and  E2  =  E^.  Suppose  I  6 
C^(J5.R),  satisfies  the  Palais-Smale  condition,  and 

(11) I(u)  =  ^{Lu,u)  +  b(u),  where  Lu  =  L\PiU  +  L2P2U  and  Li  :  Ei  — >  Ei  is 
bounded  and  selfadjoint,  where  Pi  :  E  — y  Ei  is  the  orthogonal  projection,  i=l,2, 

(12)  b'  is  compact,  and 

(J3)  there  exist  a  subspace  E  C  E  and  sets  S  C  E,Q  C  E  and  constants  a  >  uj 
such  that 

(i)  S  C  El  and  J|s  >  <1, 

(ii)  Q  is  bounded  and  I\3q  <  u, 

(Hi)  S  and  dQ  link. 

Then  I  possesses  a  critical  value  c>  a. 


Proof  of  Theorem  I. 

The  proof  follows  the  lines  of  the  proof  of  Theorem  1.1  in  [2]  by  Benci  and  Fortunate. 
We  will  first  consider  a  globally  defined  Hamiltonian  following  a  trick  introduced  by 
Rabinowitz. 

Let  t/>  G  C°°(R,R)  such  that  Jp'(s)  <  0  for  s  €  (r/3,  r/2)  and 


if  s  <  r/3 
if  s  >  r/2 


where  r  was  introduced  in  hypothesis  (H4).  We  define  the  globalized  Hamiltonian,  by 
extending  the  higher  order  term 

=  V-d  z  I  Wz.O  +  B(1  -  HI  z  D)  I  z  I’  .  (3.57) 


By  choosing  R>1  large  enough  and  using  hypothesis  (H4)  we  can  show  that 


(H,(z,  t),  z)  >  qH(z,  t)  V2  6  R2". 


(3.58) 


Also  from  hypothesis  (H4),  for  a  positive  constant  C5 

^(z,i)>cs(V’(|z|)|z|i’+H(l-tf(|z|))|zl»)  VzeR""  (3.59) 

and  from  the  definition  of  H  we  see  that  there  is  an  A/  >  1  so  that 


H(z,t)  =  R\z\\ 


Wz  G  R^",  I  2  |>  M. 


We  define  the  functional 


fk{z)  =  Q{z) 


(3.60) 


(3.61) 
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for  all  z  in  E.  It  is  well  known,  see  [13],  that  under  our  hypothesis  fk  is  well  defined  in 
E,  of  class  and  the  critical  points  of  fk  are  classical  fc-periodic  solutions  of 


i  =  JA(t)z  +  JHxi^z^t). 


(3.62) 


In  order  to  prove  the  existence  of  critical  points  for  fk  we  will  use  min-max  arguments. 
Specifically  we  will  use  Theorem  3.1  stated  above,  which  requires  that  the  linking  con¬ 
dition  given  in  (i),  (ii)  and  (iii)  be  satisfied. 

From  (3.58)  and  (3.59)  and  assuming  without  loss  of  generality  that  r  <  1  there 
are  positive  constants  cg,  c-j  and  cg  so  that 

H{z,i)>c^\z\‘l-c^,  V2GR2"  (3.63) 


H{z,t)<c,,\z\\  V2€R-". 

Now  we  can  prove  the  linking  conditions.  Let  z  E  E^ ^  then  from  (3.64) 

rife 


fk{z)=\\zfQ-  j  H{z,i)dt 
Jo 

-cs  J 


>  Z 


2  1^  dt 


=  11  2  Kq  -C8  \z{l  dt 


(3.64) 


(3.65) 


where  |  •  denotes  the  usual  norm  in  the  space  L®([0,  Ar],  R^”)  that  we  will  write  as  Ll. 
As  is  well  known,  see  [13],  for  every  s  >  1  there  is  a  constant  q(s)  so  that 


I  V  |3<  0(5)  II  u  Ilf;  .  (3.66) 

Then  (3.65),  (3.66)  and  Lemma  3.5  imply  the  existence  of  a  constant  cg  such  that 

fk{z)  <11  z  !|q  -C9  II  2  11^  .  (3.67) 

Now  (3.67)  shows  there  are  constants  p  >  0  and  a  >  0  such  that 

fk{z)>a>0,  >/zeE+,  \\z:\\q=P.  (3.6S) 


On  the  other  hand,  let  Ri  and  i?2  be  two  positive  numbers  to  be  determined  later. 
Let  (j)k  be  one  eigenfunction  of  (3.41)  associated  to  the  eingenvalue  \k  we  found  in 
Proposition  3.1.  Assume  ||  <f)k  {{2  =  1-  Let 

Qk  =  {u  +  s<f)k/  V  e  E°  +  E~,  11  V  j|Q<  R2  and  0  <  s  <  }. 
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Then  for  z  =  v°  +  v  +  s<f>k  we  have 


fk{z)=^Q(z)-  f  H{z,t)dt 

JQ 

1 

=  -  \\q +2^‘^^k  -  J  H{v  +  s(f)k,t)dt. 

If  s  =  0,  from  (3.69)  we  have 


(3.69) 


fk(z)  <  0. 


(3.70) 


On  the  remaining  part  of  dQk  we  also  have  this  inequality.  Using  (3.63)  and  Holder 
inequality  in  (3.69)  we  have 

fk(v  +  S(f)k)  <  -  II  w"  IIq  I  V  +  s<^k  I2  +C7k.  (3.71) 

By  Lemma  3.3,  and  the  orthogonality  of  E'^,E~  and  with  respect  to  <,  >2 

1 


U  +  I2  >  -  II  U  +  S(l>k  II2 
a 

>  i(ii  ■>“  111 

a 

>  i  I  II 

a 


(3.72) 


Thus,  for  a  constant  cjO  >  0 


V  +  s4>k  I2  >  Ci0(!  II  +s’). 


(3.73) 


Using  (3.71)  in  (3.73)  we  get 

/fc(w  +  s^k)  S.  2^^^k~  II  V  IIq  — cii(|  II  +.s^)  +  cjk  (3-74) 

From  (3.74)  we  find  i?i  and  R2  so  that  with  (3.70)  we  obtain 

/jt(-2)  <  0  '^zedQk-  (3.75) 

See  details  in  [13].  Now  (3.68)  and  (3.75)  imply  the  linking  condition  is  satisfied.  The 
functional  fk  also  satisfies  the  Palais  Smale  condition  as  can  be  shown  by  standard 
calculations.  Thus  fk  possesses  a  critical  point  Zk-,  that  according  to  the  min-max 
characterization  satisfies 

0<a<fk(zk)<  Coo  —  Sup{ fk{z)/  z  eQk} 
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Next  we  are  interested  in  knowing  the  behaviour  of  the  critical  value  fk{zk)  when  k  is 
large.  We  will  show  that 

0  <  fk{zk)  <  (3.76) 

with  7  =  Min{2/(p  -2),{q-p  +  2)/{q  -  2)}, 

and  C  independent  of  k.  We  note  that  because  q  +  2>p>q>  2,  the  exponent  7  is 
positive. 

Let  z  =  V  +  s(f>k,v  E  E~  0  E^,s  £  [0,i2i],  then  by  (3.59) 


fk(z)  <  -  j  H{z,t)dt 

<  V>(1  z\)\z\^  +R{1  -  ^(1  z  1)  I  H  Wi  (3.77) 

<  ~  /  V’d  2  I)  U  r  +(1  -'>P{\z\)\  Z  1«  dt. 

Let  ^  =  {t  6  [0,  A:]/  |  z(t)  |<  1}  and  B  =  {t  £  [0,  fe]/  [  z(t)  |>  1}.  Then  since  r  <  1  and 
p>q 

[  V>(|  2  1)  I  2  r  +(1  -  V’d  2  D)  1  2  1''  dt  >  /  Xa\z\p  +XS  I  2  1^  dt  (3.78) 
Jo  Jo 

where  and  xe  are  the  characteristic  functions  of  A  and  B  respectively.  Using  the 
Holder  inequality  and  p  >  q  we  have 

/  Xa\z\p  +X6  I  2  I’  dt  > 
do 

>  ;fc{2-p/2}(  f’'  I  |2  ^^)p/2  {2-9/2}  (  f'  |  p  Ji)9/2 

Jo  Jo 

>  /*  I  P  dtyl^  +  (  /" 

Jo  Jo 


(3.79) 


XBZ  P  dt)’/2). 


Case  1.  Jq  xa  I  2  P  dt  >  XB  I  2  p  dt.  In  this  situation  we  clearly  have 


We  will  distinguish  two  cases,  in  order  to  obtain  (3.76). 

his  sit 

ib{2-p/2} 

2p/2 


/  X>4  I  2  P  +XB  1  2  P  dt> 
Jo 


\2  ■ 


(3.80) 


Then  from  (3.77),  (3.78),  (3.79)  and  (3.80)  we  obtain 


fk{z)<hkS^-Cuk^^-^^^U^\^2- 


(3.81) 
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using  (3.73),  and  recalling  that  ciO  does  not  depend  on  k,  we  have 

fk{z)<Xks‘^-ci3k^^-P/^hP. 

Finally  since  Xk  =  0{l/k)  we  conclude  that 

fk{z)  <  Ci3(ci4—  —  )•  (3.82) 

Maximizing  the  right  hand  side  of  (3.82),  we  obtain  the  existence  of  a  constant  cis  such 
that 

£  1^-  (3-83) 

Czise  2.  Jq  xa  1  P  dt  <  Jg  Xb  \  z  dt.  By  the  same  kind  of  analysis  eis  Case  1,  we 
obtain  that  there  exist  a  constant  Cig  such  that 

h(^)  < 

Putting  together  (3.83)  and  (3.84)  we  obtain  (3.76). 

The  final  step  in  the  proof  is  to  show  that  the  sequence  of  solutions  we  have  con¬ 
verges  to  zero.  We  show  that  — >  0  uniformly  in  C^(R,  R^”).  This  will  finish  the  proof 

because  for  |  ^  |<  r/3  the  higher  order  terms  H  and  H  coincide.  To  begin  since  Zk 

solves  (3.62)  we  see 

-QM  +  ^j  (ffz(zk,t),Zk)dt  =  0.  (3.85) 

Then,  from  the  definition  of  fk,  (3.85),  and  using  (3.58)  and  (3.59)  there  is  a  constant 
Cl 7  such  that 

fk{zk)  >  (i  -  -)  /  (Hz{zk,t),zk)dt 
2  q  Jo 

>qd-h[  H{zk,t))dt  (3.86) 

2  q  Jo 

>  Ci7  f  V’CI  Zk  I)  I  Zk  +R{1  -  V’d  Zk  D)  I  Zk  1’  dt. 

Jo 

Thus  from  (3.76)  and  (3.86)  we  see  that 

lim  [  ip{\  Zk  I)  I  Zk  1^  +i?(l  -  rp{\  Zk  |))  |  Zk  1^  dt  =  0.  (3.87) 

k-*oojQ 

Let  us  show  next  that  ]  Zk  |oo  is  bounded,  where  |  •  |oo  denotes  the  norm  in  £°°(R,  R^”). 
Let  Ik  be  an  interval  contained  in  [0,  A:]  such  that  for  any  t  G  Ik  =  [ti  ,  72] 

lzk(t)l>M,  and  lzk(7i)l=M.  (3.88) 
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For  those  k  such  that  no  such  an  interval  exists  we  have  \  Zk  |oo<  M.  For  any  k  large 
enough  and  such  that  J*  is  not  empty  we  have  by  (3.87)  that  [0,  k]\Ik  (f>.  Since  Zk 
solves  (3.62),  for  any  t  €  A  we  have 


Zk(t)  I  <1  Zkili)  \  +  f  I  Zk{s)  I  ds 

<M+f  \  A{s)  \\ Zk(s)  \  ds  +  f  \H^(zk,s)\ds 

J  7i 

<  M  +  cu(  f  1  Zk(s)  I  +  !  Hz{zk.,  s)  I  ds 


(3.89) 


Using  (3.60)  and  (3.89) 

I  zk(t)  I  <  M  +  ci8  /  I  2A:(5)  ds 

Jyi 

<  M  +  Cig  f  Ip{\  Zk  I)  1  Zk  r  +H(1  -  tp(\  Zk  I)  I  Zk  1^  ds 
Jo 

Then  from  (3.87)  it  follows  that  |  Zk  |oo  is  bounded.  Since  zjt  solves  (3.62)  we  conclude 
that  I  Zk  loo  is  bounded.  We  finally  prove  that  this  together  with  (3.87)  implies  that 
actually  |  Zk  loo-^  0.  Let  0  <  s  <  r/3,  and  mfc(s)  denotes  the  measure  of  the  set 
{<  €  [0,  A:]/  1  Zfc(t)  |>  s},  then  from  (3.86)  we  have 


fkizk)  >  C,7  ./.(I  Zt  I)  I  Zt  I-  +JJ(1  -  D)  I  Zt  I*  dt 

>  c„Minif|..{^(i  $  I)  I  e  r  +fl(i  -  i(>(i « D)  u  I’lmnW 

>  ci7  I  5  |P  mjt(^). 

From  (3.76)  and  (3.90)  we  see  that  mk(s)  — >  0  as  A:  — >  oo.  Let  0  <  e  <  r/3,  A;  G  N  and 
t  G  R,  then  either 

I  zk(t)  1<  6  (3.91) 

or  there  is  a  G  R  so  that 

\zk{o)\^^  and  \t  —  a\<mk{e).  (3.92) 

Then,  if  K  is  such  that  |  Zk  |oo^  from  (3,92) 

I  Zk{t)  I  <1  Zk(cr)  \  +  J  I  ifc  I 

<  €  +  /\mt(e), 
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from  where  limsup;i._go  |  zjt  loo<  ^  and  since  e  is  arbitrary  we  have  |  Zk  |oo— +  0,  and 
then  since  Zk  solves  (3.62)  also  |  |oo— ♦  0  as  A:  oo.O 

Remark  3.4.  If  the  linear  system  has  one  Floquet  exponent  of  indefinite  type,  or  it 
has  one  Floquet  exponent  of  positive  type  and  one  of  negative  type  then  the  system 

z  =  JA(t)z  +  JH^(z^t) 

admits  subharmonics  in  case  H  satisfies  (H4)  or  (H4-).  This  is  not  true  in  general  for 
system  with  only  positive  or  only  negative  Floquet  exponents  as  we  will  se  in  an  example 
in  the  next  Section. 


§3.3  Benci-Fortunato’s  condition  and  a  counterexample. 

In  this  section  we  will  present  a  counterexample  to  show  that  hypothesis  (H3)  is 
necessary.  We  also  show  that  (H3’),  the  related  hypothesis  assumed  in  [2],  is  stronger 
than  ours.  In  the  process  we  obtain  a  result  on  perturbation  of  definite  matrices  (positive 
or  negative). 

Let  us  start  with  the  counterexample.  Let  D  be  the  negative  block  of  dimension 
4x4  given  by 


(000  -u\ 

0  0  -1/  -1 

1  1/  0  0  I  ’ 

1/  0  0  0  / 


(3.93) 


where  i/  >  0; 
Hamiltonian 


we  note  that  D  is  a  block  of  C3  type  with  cr  =  —  1. 
H(z,t)  =  f{t)izj  +4+4+4) 


Consider  the 
(3.94) 


with  f(t)  >  a  >  0,  a  continuous  l-periodic  function  and  z  =  (-^i ,  •2^2  5  ^3  >  •2^4)-  us 
assume  {2*}  is  a  sequence  of  solutions  of  the  system 


i  =  Dz  +  JHziz,t), 


(3.95) 


z(0)  =  z(k)  (3.96) 

and  2*  — +  0  uniformly  in  C^(R,R^”)  as  k  00.  We  will  show  that  2*  =  0,  Vk  >  N 
where  JV  is  a  certain  number  in  N.  Since  H  certainly  satisfies  (Hi),  (H2),  and  (H4), 
this  example  shows  that  (H3)  is  necessary. 

Expressing  equation  (3.95)  in  component  form  yields: 


=  -1/24  + 


dH 

dzz 


(3.97) 


22  =  -VZz  -  ^4  + 


dH 

dz\ 


(3.98) 
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h  =  tyZ2  +  Zl  - 


dH 

dz\ 


Zi  =  VZi 


dH 

dz2 


(3.99) 

(3.100) 


Let  z  =  Zit  be  a  solution  of  (3.95)  and  (3.96).  Multiplying  (3.97)  by  zs  and  (3.99)  by 
Zl ,  integrating  between  0  and  k,  and  using  (3.96)  we  obtain 


rk  fk  Qjj 

/  {vzaZz  - -^Zz)dt  -  /  (1/Z2Z1  +  Zl  -  — Zi)dt.  (3.101) 

Jq  0Z3  Jo 

Multiplying  (3.98)  by  Z4,  (3.100)  by  Z2,  integrating  between  0  and  k  and  using  (3.96) 
we  obtain 

rk  dH  dH 

/  (1/Z3Z4  +  zl  -  —Z4)dt  =  /  (1/Z1Z2  -  —Z2)dt.  (3.102) 

Jo  Jo  <^-2^2 

Subtracting  (3.101)  from  (3.102) 

dH 


r  2  dH  2  a/f  r  dH  ,  air  . 

/  2^4  -  ^^4  +  2^1  -  ^-2^1  dt  =  -  — Z2  + 

Jo  az4  azi  Jo  az2  az3 

by  (3.94)  this  implies 

f  zl  +  zldt<  /  Af{t){z\-Jr  z\)dt  <c  !\zl-\- z\f  dt 
Jo  Jo  Jo 

where  c  =  4Sup/(t).  If  zi  ^  0  or  Z4  ^  0  then  there  exist  r  S  [0,  fc]  so  that 


\  zUr)  +  zlir)  \>  -. 


Since  c  is  independent  of  k  we  obtain  that  zi  =  0  and  Z4  =  0  for  k  >  N,  where  N  is 
such  that  ]  z*(t)  |<  l/y/c  Wk>  N  Vt  €  R.  Since  zi  and  Z4  are  zero,  equations  (3.97) 
and  (3.100)  imply  that  also  Z2  and  Z3  are  zero.  Thus  z*'  =  0  for  all  k  >  N. 

Remark  3.5  In  [2]  a  counterexample  is  given  to  show  that  (H3’),  the  hypothesis  related 
to  (H3),  is  necessary.  As  we  show  later  this  condition  is  stronger  than  ours,  so  that  the 
same  counterexample  shows  the  necessity  of  (H3).  However  the  counterexample  given 
here  is  finer  because  the  linear  term  has  all  Floquet  exponents  on  the  imaginary  axis. 

Let  us  consider  now  a  block  D  of  the  canonical  form  of  a  Hamiltonian  matrix  C. 
Then  we  have  the  following  result. 

Lemma  3.6. 
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I 


If  any  Hamiltonian  perturbation  of  D  has  at  least  one  eigenvalue  on  the  imaginary 
axis,  then  D  is  of  indefinite  type. 

Proof.  Consider  the  eigenvalue  problem 

X  =  Dx  +  XJx,  (3.103) 

x(Q)  =  x(k),  (3.104) 

Since  D  satisfies  the  hypothesis  of  Lemma  2.1  in  [2],  there  is  an  N  so  that  VA;  >  AT, 
(3.103)  and  (3.104)  has  a  solution  {x^,\k)  with 

J  <  y.  (3.105) 


Let  V  =  x*(0).  Then  by  (3.104), 


V  =  exp((D  +  XkJ)k)v.  (3.106) 

The  block  D  certainly  can  not  be  of  Cl  or  C2  type  since  in  this  case  the  eigenvalue 
corresponding  to  the  block  is  outside  the  imaginary  axis,  a  small  perturbation  of  D  will 
have  all  its  eigenvalues  outside  the  imaginary  axis  too.  Let  us  assume  that  D  is  of  C3 
type.  Then  we  can  write  the  matrix  D  as 

D  =  -al3V  -  aW 

If  we  choose  a  G  [0,27r]  so  that  exp(— ia)  =  exp(—ia^k),  then  from  (3.106)  and  the 
properties  of  V  we  have 


V  =  exp{—aV  —  akW  +  \kkJ)v.  (3.107) 

By  Lemma  A. 4,  taking  the  logarithm  of  the  matrix  in  (3.107)  gives 

{—{a  +  mk)V  —  (TkW  +  XkkJ)v  =  0  (3.108) 

where  rrik  G  Z.  Dividing  (3.108)  by  k,  we  obtain  that  Xk  is  an  eigenvalue  of  the  problem 

-  aW  +  XkJ)v  =  0.  (3.109) 

Now  (3.109)  has  the  form  of  the  eigenvalue  problem  (2.2): 

—{eV-\-aW)v-\-XJv  =  Q.  (3.110) 

As  we  mentioned  in  section  §2,  the  eigenvalues  of  (3.110)  are  analytic  near  e  =  0.  Thus 
they  have  the  form 

Aj(e)  =  Aj  +  Xj€  +  AjC^  -f-  —  (3.111) 
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Noting  that  the  matrix  {el  +  oV~^W)  has  only  large  eigenvalues  for  e  large,  (3.109) 
implies  Hk  =  {—a  —  Tnk^/k  is  bounded.  Suppose  that  Hk  is  bounded  away  from  zero. 
Since  el  +  aV~^W  is  invertible  Ve  5?^  0  ,  for  e  >  cq  its  eigenvalues  will  be  bounded  away 
from  zero  contradicting  (3.105).  By  repeating  this  argument  for  any  subsequence  of 
{Hk}  we  can  conclude  that  /xfc— +0asfc-^0.  Then  for  some  j,  a  subsequence  of  {/^k}, 
and  k  large  enough 

=  A,-(/iit)  (3.112) 

If  Xj  is  the  first  nonzero  coefficient  in  (3.111),  then  (3.105)  implies 

and  then  the  block  D  can  not  be  negative. 

If  we  look  at  the  proof  of  Lemma  2.1  in  [2]  we  realize  that  there  is  an  eigenvalue 
Uk  of  (3.103)  and  (3.104)  such  that 


3tT  ^  TT 

~T  - -  ~k‘ 


Then  we  can  repeat  the  argument  given  above  to  obtain  that  D  is  not  positive.  This 
finishes  the  proof  for  a  block  of  C3  type.  For  a  block  of  C5  type  we  would  proceed  in 
the  same  way.  We  would  actually  prove  that  this  implies  a  contradiction  because  a  C5 
block  is  either  positive  or  negative.  0 

Remark  3.6  As  a  consequence  if  D  is  a  positive  (or  negative)  block  then  there  are 
arbitrarily  small  Hamiltonian  perturbations  of  D  not  having  any  eigenvalue  on  the 
imaginary  axis.  This  is  an  interesting  result  we  have  proved  in  a  very  indirect  way. 

We  now  compare  hypothesis  (H3’)  of  Benci  and  Fortunato  with  our  hypothesis 
(H3). 

(H3’)  Every  Hamiltonian  perturbation  of  C  has  at  least  one  eigenvalue  on  the 
imaginary  axis. 

Proposition  3.2 

If  a  Hamiltonian  matrix  C  satisfies  (H3’)  then  C  has  at  least  one  purely  imaginary 
eigenvalue  that  is  not  of  positive  type,  and  at  least  one  that  is  not  of  negative  type. 

Proof.  We  prove  it  by  contradiction.  If  C  does  not  have  eigenvalues  on  the  imaginary 
axis  then  clearly  (H3’)  is  false.  Let  assume  now  that  all  eigenvalues  of  C  on  the  imaginary 
axis  are  of  positive  type.  Then  by  definition,  all  blocks  of  C  either  correspond  to 
eigenvalues  outside  the  imaginary  axis  or  they  are  positive.  Lemma  3.6  applied  to  each 
of  the  latter  gives  a  perturbation  that  takes  all  eigenvalues  outside  the  imaginary  axis, 
i.e.  (H3’)  is  false.  We  could  do  the  same  argument  if  all  eigenvalues  on  the  imaginary 
axis  are  of  negative  type.  0 

Remark  3.7.  This  proposition  shows  that  our  hypothesis  (H3)  is  weaker  than  (H3’). 
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§4.  REDUCTION  TO  THE  CENTER  MANIFOLD  AND 
PROOF  OF  THEOREM  II. 

In  Section  §3  we  proved  Theorem  I  which  guarantees  the  existence  of  nontrivial 
subharmonics  converging  to  zero  for  the  Hamiltonian  system  (HS),  under  hypotheses 
(HO),  (HI),  (H2),  (H3)  and  (H4). 

In  this  section  we  will  prove  a  similar  theorem  but  we  will  impose  the  superquadratic 
condition  (H4)  for  H  only  on  the  center  manifold. 

In  order  to  obtain  our  result  we  will  need  some  extra  differentiability  on  the  Ha¬ 
miltonian  H: 

(H5)  H  is  of  class  near  2  =  0. 

Under  hypotheses  (Hi),  (H2)  and  (H5)  we  will  reduce  system  (HS)  to  a  Hamiltonian 
system 

w  =  JB{t)w  +  JK^{'w,t)  (4.1) 

on  the  center  manifold. 

We  will  replace  (H4)  by; 

(H6)  There  are  constants  c  >  0,r  >  0,q  +  2  >  p  >  q  >  2  such  that 

{Ku,iw,t),w)  >  qK{w,t)  >  c  1  u;  I’  (4.2) 

I  |<  r,  <  e  R. 

Here  m  is  the  dimension  of  the  center  manifold.  A  direct  application  of  Theorem  I  leads 
to 

Theorem  II. 

If  (H5),  (Hi),  (H2),  (H3)  and  (H6)  hold,  then  there  exists  an  integer  N  so  that  for 
every  k  >  N,  there  exists  a  nontrivial  solution  Zk  of  (HS)  with  period  kT  and 
converges  to  zero  as  k  oo  uniformly  in  C*(R,  R^”). 

The  reduction  to  the  center  manifold  allows  us  to  prove  the  theorem  under  condition 
(H6)  that  in  principle  is  weaker  than  (H4). 

As  we  will  see  later  in  this  section  the  reduced  Hamiltonian  K  is  not  known  ex¬ 
plicitly,  so  that  it  is  interesting  to  have  a  condition  on  the  original  Hamiltonian  H  that 
implies  (H6).  If  we  assume  some  hypothesis  on  the  behaviour  of  the  second  derivatives 
of  near  0  then  we  can  prove  that  (H4)  imply  (H6).  See  Proposition  4.3.  Actually  we 
only  need  to  impose  (H4)  on  the  center  manifold.  We  will  be  more  precise  about  this 
point  in  §4.4  where  the  adequate  notation  will  be  available. 

§4.1.  Reduction  to  a  system  with  a  constant  linear  term. 

With  hypotheses  (H5),  (Hi)  and  (H2)  we  can  write  the  Hamiltonian  system  (0.1) 
as: 

z  J A{t)z J (4.3) 
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where  H{z,  t)  =  o(|  z  p).  System  (4.3)  is  refered  as  (HS).  We  have  the  analog  of  Lemma 
3.1.  Using  the  same  notation  considered  there,  we  have 

Lemma  4.1. 

There  are  matrices  Cr,  R  and  P{i)  so  that  Cr  is  Hamiltonian,  R?  =  I,  and  P{t)  is 
symplectic  satisfying  P{t  +  1)  =  P(t)R,  for  which  we  have: 
z(t)  is  a  k-periodic  solution  of 

z  =  J  A{t)z  +  JHz{z,t) 

if  and  only  if  ^(t)  =  P(t)~^z(t)  is  an  Rk-periodic  solution  of 

i  =  CrC  +  m^(^,t)  (4.4) 

where 

Hi(U)  =  H{P{t)z,t).  (4.5) 

Proof.  Can  take  Cr  =  Cr,  R  =  R  and  P(t)  =  P{t)  as  in  Lemma  3.1,  and  do  the  proof 
similarly.  <C> 

Remark  4.1.  The  new  Hamiltonian  is  not  periodic  in  general,  but  it  satisfies 

+  (4.6) 

Vt  G  R,  €  R^".  We  will  say  that  Hi  is  R-periodic. 

Remark  4.2.  In  Lemma  4.1  the  matrices  Cr,  R  and  P  are  not  unique,  they  can 
be  transformed  by  symplectic  transformations.  They  can  be  chosen  so  to  fit  with  the 
setting  of  the  Center  Manifold  Theorem  presented  in  Appendix  C. 

By  definition  of  Cr  and  C,-  we  have  a  symplectic  transformation  Q  so  that 


Q-^CrQ  =  S 


and 


Q-^CiQ  =  nl 


see  (2.33)  and  before  for  the  definition  of  S  and  I.  The  matrix  S  is  a  block  matrix 


5  = 


0  S,, 


where  Sp  and  Sg  are  Hamiltonian  matrices.  Let  Qp  and  Qg  be  symplectic  matrices 
transforming  Sp  and  Sg  to  their  corresponding  canonical  forms  (See  Theorem  1.1). 
Then  if 

Q  f  Qp  0 

\0Q, 
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and 


Q2  =  QQi 


the  matrix 

Cr  =  Q2^CrQ2  (4.7) 

is  the  canonical  form  of  Cr,  and  due  to  the  block  structure  of  S  and  I 


TrI  =  Q2^CiQ2.  (4.8) 

We  define 

R  =  (4.9) 

and 

Pit)  =  Q-^Pit)Q2.  (4.10) 

For  the  matrices  Cr,  R  and  P(t)  Lemma  4.1  holds.  By  the  block  structure  of  the 
canonical  form  matrix,  we  can  reorder  Cr  so  that 


Cr  = 


0 

C2 


and 


with  Cl  having  eigenvalues  only  on  the  imaginary  axis,  C2  having  eigenvalues  only 
outside  the  imaginary  axis,  and 


Moreover  C2  and  R2  have  the  block  structure 


C2 


f  C21  0  \ 

V  0  C22  J 


and 


0 

R22 


C21  has  eigenvalues  with  positive  real  part,  C22  has  eigenvalues  with  negative  real  part 
and  the  matrices  R21  and  R22  have  the  form 
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for  t  =  2, 3.  The  following  commutativity  property  is  consequence  of  the  block  structure 
and  the  special  form  of  R: 


(^2\R2\  =  RllC^Xt  ^22^22  =  ^22^*22)  S'Hd  C\Ri  =  RiC\.  (4-11) 
Remark  4.1  holds  with  this  new  matrices. 

Remark  4.3.  In  the  canonical  form  we  considered  before  we  can  take  the  ‘off  diagonal’ 
elements  to  be  any  number  7^0.  See  [4]. 

We  can  define  symmetric  matrices  Ai  and  .42  so  that 


JjAj  —  Oj  j  —  1,2. 

Then,  after  the  transformation  we  can  write  (4.3)  as 

(4.12) 

(4.13) 

(2  =  d2Mi2  + 

(4.14) 

In  Section  §4.2  we  drop  the  index  1  in  H  and  the  bar  in  R. 


§4.2  Reduction  of  (4.13)  and  (4.14)  to  the  center  manifold. 

We  want  to  reduce  (4.13)  and  (4.14)  to  the  center  manifold.  Because  this  system 
is  not  periodic  but  R-periodic  we  will  need  a  special  version  of  the  Center  Manifold 
Theorem.  In  Appendix  C  we  provide  the  changes  necessary  to  modify  the  original 
existence  proof  in  order  to  suit  our  needs.  Our  basic  references  for  this  theorem  are  the 
original  paper  by  Kelley  [9]  and  the  book  by  Caxr  [5]. 

Theorem  4.2.  (Center  Manifold  Theorem) 

If  H  satisfies  (H5)  then  there  exists  a  function 

hiNxR—^ 

where  N  is  a  neighborhood  of  the  origin  in  R^">,  such  that 

/i(Ri6,t  +  l)  =  R2ft(6,0  (4.15) 

h{0,t)  =  0,  fi'(0,t)  =  0  (4.16) 

and  the  manifold 

is  an  invariant  manifold  for  (4.13)  and  (4.14).  M  is  called  the  center  manifold. 

Proof.  See  Appendix  C.<C> 
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System  (4.13)  and  (4.14)  can  then  be  reduced  to 


(4.17) 

or  equivalently 

(4.18) 

where  H  is  given  by 

+  2(^2<^2,if2)  +  -^(6>^2,0- 

System  (4.18)  is  a  Hamiltonian  system  if  it  is  interpreted  properly;  namely  we  have  to 
consider  in  R^”‘  a  nonconstant  symplectic  form.  In  order  to  apply  the  results  of  Section 
§3,  particularly  Theorem  I,  we  need  to  transform  (4.18)  to  put  it  as  a  Hamiltonian  system 
in  with  the  standard  symplectic  form.  Let  us  define  the  reduced  Hamiltonian 

H2Uut)  =  HUuh{^ut),t).  (4.19) 

Using  (4.6),  (4.11)  and  (4.15)  we  see  that 

■^^2(6,*  +  1)  =  -(Ai^i,^i)  +  -{A2h{(i,t  +  1),  +  1)) 

+  -^(<^1 1  <^(6  >  ^  +  1)1  ^  +  1) 

=  \{A,R,(uRi(i)  +  l{A2h{R,(ut),h{R,(,,t)) 

=  H2{Ri^ut)- 

in  other  words  we  have  proved  that 


H2Ui,t  +  l)  =  H2{Ri^ut)  VtER,  eiGiV.  (4.20) 

The  derivative  of  H2  is  given  by 

H2^,  (ei ,  t)  =  (6 ,  h{i, ,  t),  t)  +  {h'itx ,  <))"%  (6 ,  Ki, ,  t),  <)  (4.21) 

where  '  denotes  derivarive  with  respect  to  ^1. 

On  the  other  hand,  M  being  an  invariant  manifold,  the  right  hand  side  of  (4.13) 
and  (4.14)  should  be  tangent  to  M,  that  is 


(  I  =  (  A(6,0a^  +  '^'(6,0A  I  (4.22) 
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where  •  denotes  derivative  with  respect  tot,  \  E  and  /z  G  R.  From  (4.22)  we  obtain 
that 

^2%  (6 ,  ,  t),  t)  =  A'(6 ,  t)h  (fj  ,hiiu  t),  t)  +  (4.23) 

Inserting  (4.23)  in  (4.21)  we  obtain 

^2^(6,<)  =  (/  -  -  h'^J2h,  (4.24) 

and  then  system  (4.18)  can  be  written  as 

-(Ji  +  h'^J2h%  =  ^26(6,0  +  h'^J2h.  (4.25) 


This  system  is  Hamiltonian  if  we  consider  R^”*  with  the  symplectic  form  given  by 
Ji  +  /i'’’ Jo/i'.  We  will  perform  another  transformation  so  that  (4.25)  is  a  Hamiltonian 
system  in  R^"'  with  the  standard  symplectic  form. 

Before  doing  so  let  us  collect  some  properties  of  the  terms  involved  in  (4.25).  Let 
us  define 


From  (4.15)  we  see  that 
and  then 


u^i((i,t)  =  h'^((ut)J2h'Ui,t)  and 

(4.26) 

^2(6,i)  =  /i'"(ei,<)-^2A(6,<). 

(4.27) 

+  1)  =  R2h(Ri(i,t) 

(4.2S) 

+  1)  =  R2h'{Ri^\,t)Ri. 
f  i?i,  (4.26),  (4.27)  and  (4.29)  we  obtain 

(4.29) 

+  1)  =  Ri^i{Ri^i,t)Ri  and 

(4.30) 

+  1)  =  Rl<^2{Rl^l,t). 

(4.31) 

For  notational  convenience  in  the  next  section  we  will  drop  the  index  1  from 
and  ni,  and  the  index  2  from  H2- 


§4.3  Transformation  of  (4.25)  via  the  Darboux  Theorem. 

In  this  section  we  will  find  a  transformation  that  will  change  (4.25)  into  a  Hamil¬ 
tonian  system  with  constant  symplectic  form.  Before  we  go  to  the  construction  let  us 
assume  we  have 

F:V  xR  — >  W  X  R 
{w,t)  ^ »  {^,e) 

^  :y  X  R  — yR 
{w,t)  f — >  g{w,t). 
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where  F  is  a  difFeomorphism  of  class  g  is  a  function  of  class  and  V  and  V'  are 
neighborhoods  of  the  origin.  Let  us  assume  F  and  g  satisfy: 

(FO)  F(w,t)  =  (G{w,t),t) 

(FI)  G(0,t)  =  0,  G'(0,t)  =  0,  g{0,t)  =  0  and^40,t)  =  0, 

(F2)  G{w,t  +  1)  =  RG(Rw,t)  and  g(w,t  +  1)  =  g(Rw,t)^ 

(F3)  G''^{w,t){J  +  u!i(F(w,t)))G'(w,t)  =  J  and 

G''^{w,t){J  F (Ji{F{w,t)))G{w,t)  +  G'''{w,t)L02{F{w,t))  =  gw{w,t) 

where  here  and  in  what  follows  '  and  •  denote  derivatives  with  respect  to  w  and  t 
respectively. 

Lemma  4.2. 

Let  F  and  g  as  above.  Then  ^(t)  is  an  Rk-periodic  solution  of  (4.25)  if  and  only  if 


w{t),  defined  implicitly  by 

F(w(t),t)  =  m,t), 

(4.32) 

is  an  Rk-periodic  solution  of 

W  =  JH3w{w,t) 

(4.33) 

where 

H3{w,t)  =  H{F{w,t))  +g{w,t). 

(4.34) 

Proof.  Since  ^(t)  is  jRA:-periodic  it  satisfies  ((t  4-  k)  =  R^((t).  Then  from  (F2) 

R>‘({t)  =  G{w{t  +  k),t  +  k)  =  R^GiR'^wit  +  k),t) 

so  that 

^(t)  =  G{R'‘w(t  +  k),  0,  (4.35) 

and  this  implies  w(t  k)  =  R^w(t). 

Let  us  assume  now  that  ^(t)  is  a  solution  of  (4.25).  Then  differentiating  (4.32)  and 
multiplying  by  —  (7  +  wi)  we  obtain 


— (J  +u;i)^  =  —(J  +  i^i)G'(w,t)w  —  (J  +  u;i)G(u;,<).  (4.36) 

From  (4.34) 

H3wiw,t)  =  G'''H({^,t)  +gw{'w,t).  (4.37) 

Multiplying  (4.36)  by  G'’",  using  (4.25)  and  (4.37)  we  obtain 

— G  (^J  +  uji )G' w  =  H3tu{w, t)  4”  G^’'u?2  +  G'^ [J  +  u;i)G  —  gwi'^^^t),  (4.38) 


From  (F3)  we  see  that  (4.38)  is  exactly 


-Jw  =  H3y,(w,t) 
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completing  the  proof.  0 

Remark  4.4.  From  (4.20),  (F2)  and  (4.34)  we  see  that 


H3{w,t +  1)  =  H3{Rw,t),  (4.39) 

and  from  (FI)  we  have  fi'3(0,t)  =  0  and  .ff3u,(0,<)  =  0. 

It  remains  to  prove  the  existence  of  the  diffeomorphism  F  and  the  function  g.  This 
is  done  via  the  Darboux  Theorem.  We  present  here  an  elementary  proof  of  the  Darboux 
Theorem.  For  a  proof  using  the  language  of  differential  forms  see,  for  example,  [1]. 

We  will  introduce  an  extra  variable  and  find  a  function 

J'lVxRxR  — ^y'xRxR 
(uj,t,u)  I — »  F(w,t,u)  = 

Taking  appropriate  projections  of  F  we  will  define  F  and  g.  We  call  z  =  {w,t,u)  and 
Let  us  consider  the  matrices 

/J  0  0\ 

flo  =  0  0  1  I  and  (4.40) 

Vo  -1  0/ 

(  Wi(C,^)  W2(C,0)  0\ 

^^('7)=  0  0  .  (4.41) 

Vo  0  0/ 

If  we  supplement  R  with  two  extra  entries 

(R  0  0\ 

.R=  0  1  0  , 

Vo  0  1/ 

then  from  (4.30)  and  (4.31)  we  have 

fi(e,  ^  +  1,  ^)  =  m{R^.  e,  fi)R.  (4.42) 


Lemma  4.3. 


dQij 

dr]k 


d9,ki  d9,ki 
dT]i  drjj 


1  <  i,j,  k  <2n  +  2 


(4.43) 


Proof.  By  direct  calculation. 0 

Remark  4.5  Lemma  4.3  says  the  the  2-form  defined  by  Q.  is  closed. 


Because  satisfies  (4.43)  we  have  the  following  lemma. 

Lemma  4.4.  (Poincare  Lemma) 

There  is  a  function  a  :  N  x  — >  R2n+2 

*  dfy ' 

1  <  i,  J  <  2n  +  2.  (4.44) 

Moreover  the  function  a  satisfies 

a(^,e +  l,fx)  =  Ra(^,6,/j,).  (4.45) 


Proof.  We  define  a  by  the  following  integrals 

.1  2n 


<^i(^)  =  J  2^^ik(T^,9,n)T^kdT,  l<i<2n  and 


(4.46) 


fc=i 
.1  2n 


oiiiv)=  f  '^^ikiT^,6,fi)^kdT,  i  =  2n  +  l,2n  +  2.  (4-47) 

*=i 

First  let  us  assume  1  <  i,j  <  2n.  Then  differentiating  (4.46)  we  have 

2n 


5a, 

drjj 


-(ri)=^  f  flijir^,e,fx)T +  '^^^{T^,e,fi)T^^kdT  (4.48) 

>  -^0  *=i 


and 


Ip-(^)  =  /  ^ji(TLd,fi)T +  Y^^^{T(,e,fi)r^(kdT.  (4.49) 

5??,  70  jt=i 

Noting  that  =  — fiji,  subtracting  (4.49)  from  (4.48),  and  using  (4.43)  we  obtain 

~(n)  -  =  I  2n„(rf,«,^)r  +  g(^  -  ?^ymdT 


dr]j 


dr]i 


p\  <90 


(4.50) 


In  this  last  expression  we  identify  the  integrand  as 

d 


dr 


(4.51) 
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so  that 


dcxi  ,  .  dai  ,  .  ^  .  s 

g^(»)  - 


The  case  i  =  2n  +  1,  j  =  2n  +  2  is  trivial  because  Oj  =  0  ,  ai  does  not  depend  on  u  and 
Qij  =  0.  The  other  cases  are  obtained  in  a  similar  way. 

Condition  (4.45)  is  easily  obtained  by  using  (4.42)  and  the  definition  of  a.  0 

Lemma  4.5.  (Darboux  Theorem) 

There  is  a  diffeomorphism  :  V  x  >  V  x  such  that 


J^(zy^(Qo  +  n(J^(z)))J^'(z)  =  Qo. 


(4.52) 


Proof.  The  function  JF  will  be  defined  by  a  differential  equation.  For  s  6  R,  let 

^s(rf)  =  (4.53) 

Then,  since  Qq  is  invertible  and  fi(O,0,/x)  =  0,  there  is  e  >  0  so  that  ^^(j?)  is  invertible 
for  every  jj  G  ^<(0)  x  R^.  Then  we  can  define 

= -^7^(n)air}),  (4.54) 

and  consider  the  differential  equation 


T]{0)  =  2. 

Since  a{0,6,p)  =  0  y{d,p)  €  R^,  equation  (4.55)  with  the  initial  condition 

T7(0)  =  (0,<,u) 


(4.55) 


has  the  solution 

77(5)  =  (0,t,  u)  Vs  €  R. 

Then,  by  the  general  theory  of  ordinary  differential  equations  we  can  define 

H  :  B,,{0)  X  R2  X  [0,2]  — ^  R^"  x  R^ 


where  7i{z,s)  is  the  solution  of  (4.55),  and  0  <  ei  <  e.  H{z,l)  is  a  diffeomorphism  of 
class  C^,  see  [6].  We  define  ^{z)  =  '?f(2, 1).  Let  us  check  that  the  function  J-  satisfies 
(4.52).  By  definition  we  have 

0(77,  s)  =  fio  +  5^(77). 
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Then  we  have  the  following  derivatives: 


dilij 

dzk 


(n(z,s),s) 


dflij  dHi 
drji  dzk 


^  N  dQijdHi  ,  ^ 
— ^  +%• 


ds 


drji  ds 


Here,  and  in  what  follows  we  will  use  Einstein’s  convention  for  sums.  Let 
the  differential  equation  (4.55)  componentwise: 


nij(n(z,s),s)-j^(z,s)  =  -ai(n(z,s)) 

and  differentiate  it  with  respect  to  Zk-  Using  (4.56)  we  obtain 

dQ.ij  dHi  dHj  /  X  A  d  dHj  .  .  dai  dHi 

Our  purpose  is  to  show  that 

n'^(z,  i)n(n{z,  i),  i)n'iz,  i)  =  rio. 


Since  we  obviously  have 


n'^(z,o)a(n{z,o),o)n'{z,o)  =  Qo 

it  is  enough  to  show  that 

4-'H'^{z,s)n{H{z,s),s)n'{z,s)  =  0. 
as 

We  do  this  componentwise.  Differentiating  and  using  (4.57)  and  (4.58) 

d.dHi.  .  sd'Hj  , 

_r6  ^  d'Hj  dHj  ~  d  dHj  dHj  dHjdhjjdHj 

ds  dzm  dzk  ds  dzk  dzm  dzm  ds  dzk 

doLj  dHi  dHj  dai  dH,  dHi  d^ji  dHi  dHj  dHi 

—  - rL - -  -  ^  ^  - - -  ^ - £ - - 

dT]i  dzjn  dzk  drji  dzk  dzm  drji  dzm  dzk  ds 

^  dHj  dHj 

drii  dzk  dzm  ds  ^  drji  dzm  dzk  ds  dzm  dzk 


(4.56) 

(4.57) 
consider 

(4.58) 

(4.59) 


(4.60) 
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Changing  some  of  the  indices  in  (4.60)  and  reordering  the  terms  we  obtain 

drji  dr]j  dzm  dzk  ^  drji  drjj  '*'  drji  dz^  dzk  ds 

Now  (4.61)  is  zero  by  Lemmas  4.3  and  4.4. <C> 

Now  we  are  prepared  to  prove  the  existence  of  F  and  g. 

Proposition  4.1 

There  is  a  diffeomorphism  F  and  a  function  g  satisfying  (FO),  (Fl )  ,(F2)  and  (F3). 

Proof.  Let  us  consider  the  diffeomorphism  F  defined  in  Lemma  4.5,  and  let  us  write 
it  as 

F{z)  =  (G(z),  T(z),  u  +  g{z)).  (4.62) 

We  recall  that  z  =  (w,  i,  u).  From  (4.41),  (4.46)  and  (4.47)  we  see  that  X  is  independent 
of  772n+2  =  hence  G,  T  and  g  are  functions  of  (w,  t)  only.  From  (4.41)  and  (4.47)  we 
see  that  a2n+2  =  Oj  and  from  the  block  structure  of  fl,  we  have  that 


nr  = 


(J  +  swi)"^  0 

0  0 

—  SwJC*/ +  5Wl)“^  1 


0  s(J  +  SWi)  ^U}2 

0  -1 


(4.63) 


therefore  A’2n+i(^,'S)  =  0,  and  this  implies  that 


T{z)  =  t. 


(4.64) 


From  (4.42),  (4.45)  and  (4.54)  we  see  that 


x((,e  +  i,p,s)  =  RX(R{(,e,p),s). 

The  function  H{Rw,t,  u,  s)  satisfies  equation  (4.55)  with  the  initial  condition 


(4.65) 


7f(i?ic,t,u,0)  =  (i?w,<,  w), 

according  to  the  proof  of  Lemma  4.4.  On  the  other  hand,  the  function 
H»(w,t,u,s)  =  {RG{w,t  +  l,s),t,u  +  g{w,t  +  l,^),^) 


(4.66) 


(4.67) 


satisfies 


?7»(u;,t,u,0)  =  (Rw,t^u), 


(4.68) 
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and  using  (4.65) 


dT~L  d 

-^(w,t,u,s)  =  —{RG{w,t  +  +  l,n  +  g{u},t  +  l,s)) 

=  —R7{(w,  t  +  l,u,  s) 
as 

=  RX(G(w,t  +  l,s),<  +  l,u  +  g{w,t  +  l,s),s) 
=  X(RG(w,  i  +  l,s,t,u  +  g(w,  i  +  1, 3),  5) 

=  X{'Hm(w,t,u,s),s). 

Then  by  the  uniqueness  of  solutions  of  ordinary  differential  equations 

'H{Rw,t^u,  1)  =  1) 


and  hence 

G{Rw,t,u)  =  RG{w,t  +  l,u)  and  (4.69) 

g{Rw,t,u)  =  g{w,t  +  l,u).  (4.70) 

Define  F(w,  t)  =  (G(w,  t),  t).  Because  F  is  the  projection  of  X,  and  (4.64)  holds  together 
with  the  fact  given  after  (4.62),  it  follows  that  F  is  a  diffeomorphism.  It  certainly  satisfies 
(FO).  By  (4.16)  and  the  form  of  X,  F  and  g  satisfies  (Fl).  From  (4.69)  and  (4.70)  it 
satisfies  (F2).  Lastly  (F3)  is  a  consequence  of  Lemma  4.5.0 


§4.4  Analysis  of  Hypothesis  (H6). 

Our  reduction  process  led  us  to  the  Hamiltonian  system 


W  =  JHzrviw.t). 


By  Remark  4.4  we  have  that  Hj  satisfies  (Hi)  then  we  can  write  this  system  as 

w  =  J  Aiw  +  (4-71) 

with  Hz  =  o((  w  p).  A  further  transformation  of  (4.71)  in  order  to  take  it  into  a  periodic 
system  will  give 

w  =  JB{t)w  +  JKy,{w,t)^  (4-72) 

and  we  can  apply  Theorem  I  to  it  if  we  assume  (H6).  See  Section  §4.5. 

In  this  section  we  will  study  condition  (H6)  for  H3,  and  we  will  derive  some  sufficient 
conditions  for  (H6)  to  hold  in  terms  of  the  Hamiltonian  H  or  rather  Hi  from  Section 
§4.1.  For  convenience  we  will  use  the  notation  of  Section  §4.1,  and  we  will  refer  to  Hi 
as  the  original  Hamiltonian.  We  assume 
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(4.73) 


(H7)  There  are  constants  c>0,  r>0, 9  +  2>p>g>2  such  that 

(^Ux(6,o,0,6)  >  9^1(6, 0,0  >  c  Ui  r 
V6eR"‘,  Ui|<r,  vteR, 

and  a  condition  on  the  cross  derivatives 

(H8)  F,4,j,(6,0,0  =  O(|4,  r^)  and  A„4,(6.0,()  =  O(|e,  r')  near 

6=0. 

Remark  4.6.  Clearly  (H4)  implies  (H7). 

Remark  4.7.  Due  to  the  method  we  are  using  it  seems  necessary  to  have  a  condition 
like  (H8)  in  order  to  deduce  (H7)  from  (H4). 

Remark  4.8.  Under  hypothesis  (H8)  we  see  from  the  Center  Manifold  Theorem  (Ap¬ 
pendix  C)  that 

h'CiiJ)  =  0(1  6  r-")  and  A(6,<)  =  0(1  (,  I'-’).  (4.74) 

After  integrating  (4.74)  and  using  that  /i(0,<)  =  0  we  also  have 

M6,<)=0(|6  r').  (4.75) 

The  following  lemma  describes  the  behaviour  of  G  and  g  near  ly  =  0.  Let  us  write  G  as 
G{w,  t)  =  w  +  Gi(to,  t)  and  consider  ^(p)  =  Min{2p  —  4,p  —  1). 

Lemma  4.6. 

Under  hypothesis  (H8) 

G'i(u;,0  =  O(lty  (4.76) 

and 

9u,(i«,0  =  o(|  (4-77) 

Proof.  By  the  same  argument  given  in  Proposition  4.1  we  can  show  that  Ti  has  the 
form 

'H{iv,t,u.,s)  =  (G(w,t,s),t,u  +  g{w,t,s))  (4-78) 

where  naturally  G{w,t,  1)  =  G{w,t)  and  g{w,t,  1)  =  g{w,t).  To  show  the  lemma  it  will 
be  necessary  to  analyze  the  vector  field  X  given  by  (4.54).  Let  Xi  represents  the  first 
2n  components  of  X.  Recalling  that  A'  does  not  depend  on  from  (4.55)  we  see  that 
G  satisfies 

~G  =  Xi(G,t,s),  G{w,t,0)=^w.  (4.79) 

ds 

If  a  denotes  the  first  2n  components  of  o  from  (4.54)  we  have 

-^i(^.^,5)  = -(J +5W])~^a'.  (4.80) 


51 


We  have 


(4.81) 


—  (^J  ^  —  J  sJ U^i  J  (^J ')^  J  4"  ••• 

Then  from  the  definition  of  u)i  and  (4.74) 

Differentiating  (4.81)  with  respect  to  w  and  using  (4.74)  again 

(-(J  +  sa;0-*)'  =  O(|er2). 

On  the  other  hand,  frona  the  definition  of  a  in  (4.46) 

<5  =  0(1  {  p-‘) 

and  differentiating  a  with  respect  to  ^ 

a'  =  0(1  e  IW-') 

Then  from  (4.80),  (4.82)-(4.85) 

Xi  =  o(!  ( pp-3) 

and 

X[  =  0(1  {  I'X')). 

Differentiating  (4.79)  we  find  that  G'  satisfies  the  differential  equation 

^G'  =  X[  G'  and  G'{w,  t,  0)  =  I. 
as 

Then,  integrating  (4.88)  and  using  (4.87)  we  have 

\G'{w,t)-I\<\  f\[G'ds\<cf  \X[\ds<ci  f  \G{w,t,s)f^PUs. 
Jo  Jo  Jo 

By  (4.87)  we  see  that  X[{0,t,s)  =  0,  then  noting  that  (7(0, t,s)  =  0  we 
G'(0,t,s)  =  I,  from  where 

G{w,t,s)  =  lu  +  o(l  u;  )). 


Thus  from  (4.89)  and  (4.90) 


G'{w,t)-I\<C2\w\I^^P^  . 


(4.82) 

(4.83) 

(4.84) 

(4.85) 

(4.86) 

(4.87) 

(4.88) 

(4.89) 

find  that 

(4.90) 

(4.91) 
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To  show  (4.77)  we  consider  the  2n  -f  2  component  of  (4.55),  then  we  have  that  g  satisfies 


~g{w,t,s)  =  X2n+2iG{w,t,s),t,s). 
But  from  the  definition  of  X  and  Q 


(4.92) 


^2n+2  =  (-5a;2(J  +  swi)  \1,0)q;. 


Here  u}2  =  h''^Jh  so  that 


=  0(1  (  p"-’). 

Since  h'(0,t)  =  0  we  have  A'(^,t)  =  o(l)  then,  using  (4.74) 

4  =  0(1  c  r-^). 

Now,  using  the  definition  of  a,  for  the  2n  +  1  component  we  have 

02„+i  =  o(ur^) 

and  then,  using  (4.81)-(4.S5),  (4.95)  and  (4.96),  from  (4.93)  we  find  that 

.Yk+2  =  0(1  r'). 

Differentiating  (4.92),  and  doing  an  analysis  similar  to  the  one  above  we  obtain 

9w{w,t)  =  o(l  w  j'’“^).C> 


(4.93) 

(4.94) 

(4.95) 

(4.96) 

(4.97) 


Proposition  4.2. 

If  Hi  satisfies  (Ho), (Hi ),  (H2),  (H7)  and  (H8)  then  there  are  constants  f  >  0,g  > 
2,c  >  0  so  that  q  +  2>p>q>2  and 

(H'3u;(io,t),u;)  >  qHz(w,t)  >  c  |  a; 

Vu»  €  R^",  I  ta  |<  f,  Vt  G  R. 


Proof.  Let  us  consider  the  following  functions 

a{w,t)  =  (>liia,Gi(ic;,f))  +  ^{AiGi{w,t),Gi{w,i)),  (4.98) 

K^,t)  =  ~{A2h{G{w,t),t),h{G{w,t),t)),  (4.99) 

c{w,t)  =  Hi{G{w,t),h{G{w,t),t),t),  (4.100) 

53 


Then 


H3{w,t)  =  a{w,t)  +  b{w,t)  +  c{w,t)  +  g{w,t). 

We  analyze  each  term  separately.  Differentiating  (4.98)  and  taking  the  inner  product 
with  w  yields 

+  (Aiw,G'i(w,t)w)  +  (AiGi(w,t),G[(w,i)w).  (4.101) 

From  Lemma  4.6  we  have 

(a'(u;,<),u;)  =  0(\  w  +  0(\  w  (4.102) 

By  the  definition  of  /3(p),  and  noting  that  p  >  2 

I3{p)  +  2  =  Min{2p-  2,p  +  1}  >  p  +  ^ 
where  8  =  p  —  2.  Also  2/3(p)  -f  2  >  p  +  5,  so  that 

(a'(u;,i),  ly)  =  0(1  ti)  (4.103) 

Differentiating  (4.99)  and  taking  the  inner  product  with  w  we  obtain 

=  (A2h{G{w,i),t)iG''^  h'{G{w,i),t)iv). 

Then  from  Lemma  4.6 

(6'(ty,t),iy)  =  0(j  ty  1^“^'^).  (4.104) 

Finally  let’s  analyze  the  function  c.  Differentiating  (4.100),  and  taking  the  inner  product 
with  w 

{c'{w,t),w)  =  {H^^{G{w,t),h{G{wA),t)A),G'{w,t)w) 

(4.10o) 

+  (Hi^j{G(w,  t),  h(G(w,  t),  t),  t),  h'(G{w,  t),  t)G'{w,  t)w) 

By  Taylor’s  theorem  and  Lemma  4.6  and  using  that  =  o(l)  we  have 

Hi(G{w,t),h{G(w,t),t),t)  =  +  o(|  w  |P~‘).  (4.106) 

In  a  similar  fashion  we  obtain 

{Hi^^{G{w,t),h{G{w,t),t),t),h' {G{w,t),t)G' {w,t)w)  =  o(|  w  )'’).  (4.107) 

From  (4.106),  (4.107)  and  (4.104) 

{c'{w,t),w)  =  (.ffi{,(u),0,t),u;)  +  o(j  lu  j'’).  (4.10S) 


54 


From  (4.77),  (4.103),  (4.104)  and  (4.108)  we  obtain 

+  o(|  w  |P).  (4.109) 

By  analyzing  a,  b,  c  and  g  in  the  same  way  we  also  obtain 

jH(w,0,t)  =  H3{w,t)  +  o(|  w  1^).  (4.110) 

From  (H7),  and  (4.109)  we  have 

{H3u,iw,i),w)  >  qHi{w,Q,t)  +  o(|  w  j^),  (4.111) 

thus  by  (4.110),  (4.111)  and  (H7)  there  is  e  >  0  so  that  for  ri  >  0  small  enough 

(H3y,(w,t),w)  >  {q  -e)H3{w,t),  (4.112) 

Viy  G  R^",  I  w  (<  Ti,  Vt  6  R. 

The  constant  e  >  0  can  be  chosen  so  that  q  =  q  —  e  satisfies  q  +  2>p>q>2.  On  the 
other  hand  from  (H7)  and  (4.110) 

H3(w,t)  =  Hi(w,0,t)  +  o(|  10  1^)  >  c  I  ID  1^  +o(|  w  1^) 

so  that  there  is  a  constant  rz  >  0  so  that  for  a  constant  c  and  \  w  \<  r2 

(q  -  €)H3(w,t)  >  c  1  to  |P 

by  setting  f  =  Min{ri,r2}  we  complete  the  proof.O 

Remark  4.9  As  we  said  before,  all  the  results  we  have  proved  in  this  section  can  be 
obtained  if  we  assume  that  H  satisfies  an  adequate  hypothesis.  We  preferred  to  impose 
the  conditions  on  Hi  instead,for  simplicity. 

§4.5  Proof  of  theorem  II. 

This  is  a  direct  application  of  Theorem  I  noting  that  by  Remark  4.4  H3  satisfies 
(HI).  We  only  want  to  mention  that  in  the  case  the  system  has  —  1  eis  a  Floquet 
exponent,  the  reduced  system  (4.33),  or  system  (4.71)  obtained  from  (4.33)  by  applying 
Proposition  4.2: 

w  =  JAiw  +  JH3w{w-,i) 

is  i?-periodic.  A  change  of  variables  we  perform  now  allows  us  to  obtain  (4.72)  to  which 
theorem  I  is  applicable. 

Lemma  4.7. 
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There  is  a  symplectic  matrix  function  R(t)  so  that  w(t)  is  an  Rk-periodic  solution 


of 

if  and  only  if  z(t)  =  R{t)w{t) 


w  =  JAiw  +  Jilzu){w,t) 
is  a  k-periodic  solution  of 


z  =  JB(t)z  +  Jkz{z,i) 


with  the  matrix  B{t)  defined  by  B(t}  =  -JR(t)R~\t)  +  R-'^^(t)AiR-^(t)  symmetric 
and  the  higher  order  term  given  by  K(z,t)  =  H3(R{t)~^  z,t). 

Proof.  The  proof  corresponds  to  the  nonlinear  reduced  version  of  Lemma  3.2.  The 
matrix  R  we  axe  using  here  corresponds  to  the  block  Ri  we  defined  in  Section  4.1.  We 
can  define  R(i)  in  the  same  way  as  we  did  in  Lemma  3,2.  The  proof  follows  directly.O 


56 


APPENDIX  A.  LOGARITHM  OF  A  MATRIX. 


This  appendix  is  devoted  to  the  definition  of  the  logarithm  of  a  real  matrix.  We 
collect  here  some  elementary  properties  of  the  logarithm,  as  they  are  used  in  the  text. 
We  will  define  the  logarithm  of  a  matrix  using  the  Cauchy  formula  for  the  logarithm 
in  the  complex  plane.  We  will  follow  the  presentation  of  YaJcubovich  and  Starzhinskii 
[16],  and  we  refer  to  this  book  for  the  proofs  we  omit. 

Let  X  be  an  n  X  n  matrix  with  coefficients  in  C.  Let  Ai,  A2, ...,  A^  be  the  different 
eigenvalues  of  X  (without  counting  multiplicity).  We  assume  that  X  is  nonsingular,  so 
that  all  eigenvalues  axe  different  from  zero.  Let  /  be  a  ray  in  the  complex  plane  going 
from  the  origin  to  infinity,  not  containing  any  of  the  eigenvalues  Ai,  A2, ...,  A*,  and  not 
containing  1.  In  C  \  /  we  can  define  a  single  valued  branch  of  the  logarithm  in  the  usual 
way 

log  z  =  log  I  z  I  arg  2  (^-1) 

where  axg  z  is  the  branch  of  the  argument  function  satisfying  arg  1  =  0. 

Let  Fj  be  circles,  centered  at  the  points  \j,j  =  1, 2,  ...k,  which  are  mutually  disjoint 
and  do  not  intersect  /.  Let  Gj  denote  the  open  disk  with  boundary  Fj.  If  G  =  ljj_i  Gj 
then  log  2  is  well  defined  and  analytic  in  an  open  set  containing  G. 

Using  the  Cauchy  formula  we  now  define  a  logarithm  function  for  the  matrix  X  by 

losX  =  :^^'£j^(zIn-X)-'\otzdz.  (A.2) 


Lemma  A.l. 

(i)  If  K  =  logX  then  X  =  exp  K. 

(a )  If  Q  is  a  nonsingular  matrix  then 

\ogiQ-^XQ)  =  Q-\\ogX)Q. 


Proof.  For  (i)  see  [16],  page  55.  (ii)  is  direct  from  the  definition. <C> 

Remark  A.l.  A  more  general  notion  of  logarithm  can  be  defined  if  in  each  Gj  a 
different  branch  of  the  complex  logarithm  is  considered.  If  mj  G  Z,  then  define 

log(2)my  =  logz  +  2irimj  (A.3) 

then  we  define 

\og{XU  =  —  E  /  X)-^  log{z)m,  dz,  (.4.4) 

27rz  Jtj 
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where  m  =  (mi,m2,  This  logarithm  of  X  also  satisfies  Lemma  A.l.  For  our 

purposes  it  will  be  enough  to  take  rrij  =  0,_;  =  1, 2, k.  See  [16]. 

Remark  A. 2.  If  the  matrix  X  is  real,  one  would  like  logX  to  be  real.  This  is  not 
always  possible,  however  the  following  lemma  tell  us  when  it  is  so. 

Lemma  A. 2. 

If  X  is  3  real  matrix  with  no  eigenvalues  on  the  negative  real  axis  then  logX  is  a 
real  matrix. 

Proof.  See  [16],  page  56. 

When  the  matrix  A”  is  a  symplectic  real  matrix,  i.e. 

X^JX  =  J 

then  the  logarithm  of  X  satisfies  a  special  property: 

Lemma  A. 3. 

If  X  is  a  symplectic  matrix  then 

(logXyj  =  -  JlogX.  (A.5) 


Proof.  See  [16],  page  211. <> 

Remark  A. 3.  In  (A.5)  we  have  used  *  for  conjugate  transpose  instead  of  r  because  in 
general,  as  we  said  before,  the  matrix  log  A  is  not  real. 

The  block  property  for  the  logarithm  of  a  matrix  is  obvious,  but  we  make  it  explicit. 
Let  the  matrix  A  have  the  following  block  form 


A  = 


where  A,-  is  an  n,-  x  n,  matrix,  i  —  1,2.  Let  us  further  assume  that  the  eigenvalues  of 
Ai  are  different  from  the  eigenvalues  of  A2,  then  we  obtain  directly  from  the  definition 


that 


(A.6) 


Here  we  understand  the  logarithm  of  the  blocks  defined  as  in  (A. 2)  but  restrict  the  sum 
to  the  eigenvalues  in  the  corresponding  block. 


Lemma  A. 4. 

If  X  is  as  before  and  both  log(— A)  and  log(A')  are  well  dehned  then 


log(-A)  =  in  In  +  log  A 
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Proof.  For  every  z  €  C  such  that  z  and  —z  do  not  belong  to  I  we  have 


log(— 2)  =  tir  +  log  z. 


Changing  the  variable  z  to  —z  in  the  integral  (A. 2),  the  result  follows. 

Finally  we  have  a  lemma  concerning  the  relation  between  the  eigenvalues  of  a  matrix 
and  its  logarithm. 

Lemma  A. 5. 

Let  A  be  a  real  matrix.  If  exp  A  has  one  eigenvalue  equal  to  1  then  3m  €  Z  so  that 
A  —  2mV'  has  zero  as  an  eigenvalue.  Here  the  matrix  V  is  defined  by 


V  = 


and  the  matrix  K  is  defined  by 


K  = 


n 


0  IN 
1  0 


\1  ...  0  0/ 


Proof.  Let  J  be  the  Jordan  canonical  form  of  A,  i.e.  for  some  nonsingular  matrix  Q 

QJQ-^  =  A. 

Since  exp  A  has  1  as  an  eigenvalue  3v  €  R",v  ^  0  such  that 

exp  Av  =  V 

and  then 

exp(pA)v  =  u,  Vp  G  Z. 

Using  the  Jordan  canonical  form  we  obtain,  by  taking  w  =  Q~^v 

exp(pj^)w  —  w  Vp  6  Z, 

and  by  blocks 

exp{pjj)wj  =  Wj  VpG  Z,j  =  l,2,...,fc. 

Since  u  ^  0,  for  at  least  one  index  the  vector  wj  is  not  zero.  Considering  this  index, 
but  droping  it  from  the  formulae  we  have 

exp(pj)w  =  w  Vp  G  Z,  (A. 6) 
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with  J  =  \I  +  U  and  17  is  a  upper  diagonal  matrix.  Actually  the  matrix  U  has  the 
form 


where 


/U, 

0  .. 

0  \ 

0 

U2  .. 

u 

^  0 

0  .. 

UkJ 

1  0  . 

..  0\ 

0 

0  1  . 

..  0 

0 

0  0  . 

..  1 

\0 

0 

0 

..  0/ 

If  A  =  a  +  then  (A. 6)  can  be  written  as 


exp(ip/3I  +  pU)w  =  w 


Since  this  has  to  be  true  Vp  €  Z,  and  w  0,  we  conclude  that  o  =  0.  Consequently  A 
is  purely  imaginary. 

Using  the  block  structure  of  U  we  have  that  for  some  Uj,  the  corresponding  vector 
IV j  ^  0  and 

exp{pUj)wj  =  Wj,  'ip^Z  {A.Tj 

For  simplicity,  let  us  drop  the  index  j  in  (A. 7).  We  can  calculate  exp{pU)  explicitly 


(I 


exp{pU)  =  . 


P 

1 


ei 

2! 


P 


0  0  0 
\0  0  0 


^ ) 


Thus 

p” 

exp(pU)w  =  {wi  +  pW2  +  ...  ^ - rWn,  ■.■,Wn). 

nl 

All  components  of  exp{pU)w  have  to  be  bounded,  in  particular  the  first  one,  and  this  is 
possible  only  if  W2  =  W3  =  ...  =  «;„  =  0,  consequently  w  =  (u?i,  0, ..., 0),  and  then  (A. 7) 
implies 

=  1  Vp  €  Z, 

in  particular  =  1,  thus  /?  =  27rm,  for  some  m  G  Z. 

We  conclude  that  A  has  an  eigenvalue  A  =  27rim,  for  m  E  Z. 
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If  m  =  0,  then  certainly  A  —  2-KmV  has  zero  as  £in  eigenvalue.  If  m  ^  0,  since  A  is 
real,  A  =  — 27rim  is  also  an  eigenvalue  of  A.  Then  the  conclusion  follows  since  the  real 
representation  of  the  pair  (A,  A)  is  the  matrix 

/  0  —2irm  'N 

V27rm  0  J  ’ 

then  A  —  2TvmV  has  zero  as  an  eigenvalue.  The  multiplicity  of  0  is  in  this  case  at  least 
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APPENDIX  B.  BRIEF  PERTURBATION  THEORY 
FOR  SYMMETRIC  MATRICES. 


In  this  appendix  we  collect  some  basic  facts  about  the  perturbation  theory  for 
eigenvalue  problems  for  symmetric  matrices  and  we  give  a  proof  of  Propositions  2.1,  2.2 
and  2.3. 

The  subject  of  perturbation  of  eigenvalue  problems  is  well  known,  and  we  refer  the 
reader  to  the  book  by  Kato  [8]  for  more  information.  Also  the  book  by  Rellich  [14]  is 
worth  mentioning. 

After  enunciating  the  basic  theorem,  we  describe  the  calculations  necessary  to  ob¬ 
tain  the  power  series  for  the  eigenvalues  of  the  perturbed  problem.  We  then  specialize 
the  results  to  perturbation  linear  in  e  and  then  apply  these  results  to  our  situation. 

Suppose  we  have  a  power  series 

A(€)  =  Ao  -b  eAi  +  e^A2  +  ...  (5.1) 

which  converges  for  small  e,  where  the  Ai  axe  n  x  n  symmetric  matrices.  We  consider 
the  eigenvalue  problem 

A{e)v  =  Av.  (5.2) 


Theorem  B.l. 

Suppose  A  is  an  eigenvalue  for  the  matrix  A{0)  =  Aq  of  multiplicity  k  >  1,  and 
suppose  that  the  interval  (A  —  a,A  +  /3)  contains  no  eigenvalues  of  Aq  other  than  A, 
where  a,  0  are  positive  real  numbeis. 

Then  there  exist  convergent  power  series 


A,(e)  =  A  -b  A^  -b  A^  -b  ...,  i  —  1, 2, ...,  k 

(5.3) 

and 

‘>^i(e)  =  +  vp^e-b  -b  ...,  i  =  1,2, 

(5.4) 

where  A^"^^ 

6  R,  €  R”,  such  that 

A(e)v{(e)  =  Xi(e)vi(€)  i  =  1, 2, ...,  k 

(5.5) 
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and  {vi(e),vj(€))  =  Sij  =  (B.6) 

For  0  <  a'  <  a  and  0  <  <  ^0,  there  is  S  >  0  such  that  in  the  interval  (A  —  a',X  + 

Ai(e)  are  the  only  eigenvalues  of  A(e)  for  |  e  |<  ^. 

For  a  direct  proof  of  this  theorem  see  [14].  See  also  [8]. 

Our  aim  now  is  to  give  a  formula  to  compute  the  coefficients  A^"^^  of  the  power  series 
(B.3).  Given  a  symmetric  matrix  A,  we  can  consider  its  spectral  decomposition 

m 

^  =  ^AiPi  (5.7) 

«=i 

where  {Ai,  A2, ...,  Am)  =  <7(A)  is  the  spectrum  of  A  and  Pi  is  the  matrix  of  orthogonal 
projection  onto  the  eigenspace  associated  to  A,-. 

The  resolvent  of  the  matrix  A  is  defined  for  every  ^  G  C  \  a(A)  by 

R{()  =  {A-(ir'  (B.8) 


We  can  obtain  a  power  series  expansion  about  an  eigenvalue,  say  Ai,  using  (B.7):  for 
large  G  C 

^  y=o  ^ 


j=l  ^  i=l 
m  1  ^  \ 

t=i  ^  ^ 

m 


(B.9) 


i=l 


where  we  used  the  facts  that  the  series  is  absolutely  convergent  and  that  the  projections 
are  orthogonal.  We  observe  that  (B.9)  mahes  sense  for  every  ^  ^  o’(A). 

Reordering  terms,  for  ^  close  to  Ai,  we  have 


m  -  m  -  00  ^ 

1=2  *  ^  1=2  7=0  ^  ' 
00  m  - 


(B.IO) 


;=0  t=2 

Then  we  can  write,  from  (B.9)  and  (B.IO) 


(B.ll) 


j=0 
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where  the  matrix  S  is  called  the  reduced  resolvent  and  it  is  given  by 


From  (B.ll)  we  obtain  that 


(B.12) 


where  F  is  a  small  circle  around  Ai.  Considering  the  matrix  valued  function,  A(e),  we 
can  define  the  resolvent 

=  (B.13) 

If  A  is  an  eigenvalue  of  A(0)  of  multiplicity  k,  then  by  the  perturbation  theorem, 
for  small  e.  in  an  open  interval  about  A  the  only  eigenvalues  of  A(e)  are  Ai(e), ...,  Ajt(e). 
If  we  denote  by  P(e)  the  projection  onto  the  eigenspace  associated  to  Ai(e), A;t(e) 
then  using  essentially  the  argument  given  above,  we  can  prove  that 


J  R(Mdi  =  P{e)  (B.M) 

where  F  is  a  small  circle  containing  Ai(e), ...,  Ajt(e).  Noting  that 

(A(e)  -  e)  =  /  +  (^  -  e), 

from  (B.14)  we  obtain  that 

(^(6)  -  A/)P(6)  =  _  0  di.  (P.15) 


On  the  other  hand,  if  B{t)  =  A,e*, 

i2(^,e)  =  (Ao+B(e)-ar' 

=  (Ao  -a)-'(/  +  B{6){Ao  - 


Then,  a  power  series  expansion  leads  to 


p=0 
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Finally,  from  (B.15)  and  (B.16) 

{A(e)  -  A/)P(e)  =  _  X)R(()(B(X)R{()Y  d(.  (5.17) 

Knowing  the  power  series  B{e)  allows  to  obtain  a  power  series 

OO 

(X(e)  -  A)P(e)  =  Y, 

p=l 

where  the  first  integral,  for  p  =  0,  vanishes  by  (B.15).  Dividing  (B.18)  by  e  we  obtain 

1  °° 

=  -{A(€)  -  X)P(,)  =  Y 

p=0 

The  matrix  series  is  symmetric  and  for  every  eigenvalue  of  Ai  we  can  use 

the  perturbation  theorem  to  obtain  a  power  series  of  eigenvalues  of  Thus  the 

coefficients  A|^\  1  <  i  <  fc,  in  (B.3)  correspond  to  the  eigenvalues  of  Ai. 

This  process  called  Reduction  Process  can  be  continued  to  obtain  the  higher  order 
coefficients  of  (B.3).  The  matrix  .4^^^(e)  given  in  (B.19)  is  a  perturbation  of  Ai  = 
A^^^O);  thus  for  each  of  the  eigenvalues  A^^  of  Ai  we  can  apply  the  Reduction  Process 
to  obtain  the  coefficients  A^^\  By  applying  it  several  times  we  can  obtain  any  coefficient 
in  (B.3). 

So  far  we  have  described  the  general  perturbation  problem.  We  now  specialize  to 
linear  perturbations  in  c.  In  particular  we  are  interested  in  obtaining  a  formula  for  Ap 
in  the  case  all  Ai  are  zero  for  i  <  p.  This  will  give  a  formula  for  the  first  nonzero 
coefficient  in  (B.3). 

If  A(e)  =  Ao  +  cAi  then  B(e)  =  cAi  and  formula  (B.17)  greatly  simplifies: 

A,={-ir*'^J^(R{o(A,R(()r'i(,  {B.20) 

where  we  assumed,  for  simplicity  that  A  =  0.  The  resolvent  (B.ll)  is  given  by 

fl(«  =  -ip  +  f;(0>S^+'  (5.21) 

^  >=o 

If  we  rewrite  R{i)  as 

OO 

m  =  Y  (^-22) 

7=-l 
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with  the  obvious  identification  with  (B.21)  to  define  the  matrices  then  we  see  that 
given  p 

ip=  (5.23) 

where  Clp  =  {u  =  {uq,  r'l, Vp)/  Uj  >  —1,  ~  ~2}.  Thus  for  example, 

Ai  =  PA\P^  and 

A2  =  -PAiPAxS-  PAiSAiP  -  SAiPArP. 

Formula  (B.23)  gives  the  matrix  Ap  in  terms  of  P,  S  and  Ai,  for  any  positive  integer  p. 
We  are  interested  only  in  the  first  nonzero  matrix  Ap,  and  in  this  situation  the  formula 
is  even  simpler. 

Lemma  B.l 

If  Aq  =0  for  all  1  <  q  <p  then 

Ap+i  =  (p  24) 


Proof.  Our  proof  is  based  on  the  following  simple  combinatorial  fact  whose  proof  we 
present  later: 

for  5  6  N,  if  1/  €  O5  and  u  ^  (— 1,0,  ...,0,  — 1)  then  there  is  j  such  that  0  <  j  <  q 
for  which  ~ 

We  prove  the  lemma  using  an  induction  argument  on  p.  If  p  =  1,  then  Ai  =  0 
implies  A2  =  —PA^SAiP  as  we  can  easily  see  from  the  formulae  given  above.  Now 
assume  that  the  result  of  the  lemma  is  true  for  every  r  such  that  1  <  r  <  p  —  1,  and 
let  us  assume  that  Ap  =  0  for  all  1  <  g  <  p.  Let  v  €  flp+i  and  let  us  assume  that 
V  ^  (—1,0,  ...,0,  —1).  Then  from  our  claim,  there  is  0  <  j  <  p  such  that 

Qiu)  = 

=  (P^‘'«>AiP(''‘>...AiP(''^>)Ai(P(‘'>+>\..AiP(‘''+‘)) 

where  i)  ~  ~  Case  i)  We  have  two  possibilities,  a) 

(i/o,...,i/j)  =  (-1,0,  ...,0, -1)  or  b)  (t'o,...,i/j)  =  (-1,0,  ...,0, -1).  If  a)  by  induction 
hypothesis  P^‘'‘>^Ai  ^  =  Aj.  Consequently,  from  our  assumption,  (B.24)  shows 

Aj  =  0  and  hence  Q(t^)  =  0.  If  b)  then  for  =  {vo,  and  Q{v^)  = 

5(*'o)AiP(‘'i)...AiP(''i)  we  can  iterate  the  argument.  We  note  that  j  <  p  +  1.  Case  ii) 
We  proceed  as  in  Case  ii).  After  a  finite  number  of  steps  we  will  find  Q{i^)  =  0.  Formula 
(B.23)  finally  finishes  the  proof. 

Now  we  prove  the  claim  we  made  at  the  beginning.  We  have  three  cases;  i/q  =  — 1, 
1^0  =  0  or  i/Q  >  0.  Case  1.  If  i/q  =  —1  then  either  i/j  =  ...  =  =  0  which  is  not 
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allowed,  or  there  is  k  <  q  —  1  such  that  vi  =  ...  =  i^k-i  =  0  and  uk  ^  0.  If  Uk  =  — 1 
we  can  take  j  =  k.  If  i/jt  =  1  then  again  we  can  take  j  =  k.  The  remaining  situation 
corresponds  to  Vi  >  0.  Now,  since  X)i=o  ~  —  “1’  there  will  be  j  <  g 

such  that  Case  2.  If  i/q  =  0  then  j  =  0.  Case  3.  If  i/q  >  0  then  Yli=o  ^ 

for  k  =  0,  and  we  proceed  as  we  did  in  Case  1. 

In  what  follows  we  apply  the  perturbation  theory  as  presented  above  to  prove 
Propositions  2.1,  2.2. 

Proof  of  Proposition  2.1.  First  we  identify  the  terms  according  to  the  notation  used 
in  the  Appendix. 


Let 


T  =  oWJ  and  T^^'>  =  VJ. 

(P.25) 

Then  (2.3)  can  be  written  as 

(r  +  eT(^^)u>  +  A«;  =  0. 

(P.26) 

We  apply  the  reduction  process  to  (B.26). 

We  first  compute  the  matrix 

pCi)  _  pp(i)p 

(P.27) 

where  P  is  the  projection  matrix  onto  the  eigenspace  associated  to  the  eigenvalue  0  of 
T.  Displaying  the  actual  form  of  T  we  see  that  the  eigenspace  associated  to  zero  is  span 
{^1,62/}  so  that  P  takes  the  simple  form 

P  =  diag(l,0,...,0,l).  (P.28) 

From  the  definitions  given  in  (1.9),  (1-19)  and  (B.25)  we  have 

r<'>  =  ("^  “  )  ,  (S.29) 

and  then  the  coefficients 

tII’  =  T<l’,  =  0.  (B.30) 

From  (B.28)  and  (B.30)  we  see  then  that 

f(i)  ^  pr(i)p  ^  0.  (P.31) 

This  indicates,  according  to  the  reduction  process,  that  Aj  =0,  j  =  1,2,  in  (2.5). 

The  next  step  in  the  reduction  process  requires  the  introduction  of  the  reduced 
resolvent  of  T 

5  =  E  (B.32) 
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where  \j  is  an  eigenvalue  of  T  and  Pj  is  the  projection  matrix  onto  the  eigenspace 
associated  with  the  eigenvalue  Xj. 

In  order  to  compute  S  we  first  obtain  the  nonzero  eigenvalues  of  T.  After  a  calcu¬ 
lation  one  obteiins  that 

=  (B.33) 

The  eigenvalues  of  T  different  from  zero  are  a  and— cr.  Since  cr  takes  the  values  1  or  — 1, 
from  (B.32)  and  the  spectral  decomposition  of  T  we  see  that 


5  =  r.  (5.34) 

For  a  given  p  G  N  we  compute  ,  from  (B.25)  and  (B.34) 

(5j(i))P-i  =aP~i(^\vJVjy-\  (B.35) 

From  the  definition  of  V  and  W  given  in  (1-19)  we  have 

so  that,  from  (B.35)  and  (B.36) 

(57.(1)),-.  ^  ^,-1  ^  3,.j 

Then,  from  (B.29)  and  (B.37)  we  obtain 

7(l)(ST(i)),-i  =  (*■<''■’  )  .  (B.38) 


If  p  <  /  then,  using  Lemma  B.l  and  (B.38) 

fip)  =  (_i)P+ip2’(i)(5r(i))P-ip  =  0, 

and  if  p  =  /  then,  using  Lemma  B.l  and  (B.38)  again 

j>(/)  _  ^_j^(J-n)pj'(i)^52’(i)y-ip 
=  (-l)'+^<T'“Miag(l,0,...,0,l). 


(5.39) 


(5.40) 


Thus,  the  first  nonzero  coefficients  in  (2.5)  are  the  nonzero  eigenvalues  of  so 


that 


A<'>  =  (-1)V-',  ;  =  1,2. 


(B.41) 
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This  finishes  the  proof.  <C> 


Proof  of  Proposition  2.2.  In  this  case 


"o')’  ““  r<»  =  -.x 


(B.42) 


T  has  A  =  0  as  an  eigenvalue  of  double  multiplicity,  with  eigenvectors  ej  and  621- 
The  projection  matrix  is  given  by  P  =diag(l,0,  1),  and  as  before 

j(i)  _  pj-lilp  _  Q 

We  compute  the  reduced  resolvent  5.  A  calculation  gives 

det(r  -  XI)  =  A2(A2  -  l)'-i ,  (5.43) 

thus  the  eigenvalues  of  T  are  A  =  0,  A  =  1  and  A  =  —  1,  both  1  and  —1  of  multiplicity 
/-I. 

From  (B.32),  and  the  spectral  decomposition  of  T  we  see  that 


S  =  T. 


From  (B.42)  and  (B.44)  we  get 


=  i 


—d  o 
0 


Then,  for  p  6  N 


and  from  (B.42) 


(5.44) 


(5.45) 


(5.46) 


Thus  if  p  <  I, 


pCp)  _  (_i)p+i pp(i)(5p(i)^?-ip  —  Q 


and  if  p  =  / 


(d-)p-l ) 

(5.47) 

-^P  =  0, 

(5.4S) 

0,...,0,1). 

(5.49) 

Now,  recalling  that  in  this  case  I  is  odd, 

=  g(-l,0,...,0,l). 
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whose  nonzero  eigenvalues  are  1  and  —1,  giving 

=  1  and  A^^^  =  —  1. 

from  which  (2.10)  follows.  ^ 


APPENDIX  C.  THE  CENTER  MANIFOLD  THEOREM. 


In  this  appendix  we  provide  the  modifications  necessary  to  be  made  to  the  standard 
Center  Manifold  Theorem  in  order  to  suit  our  needs. 

We  consider  the  following  system  of  ordinary  differential  equations 


6  =  \ 

X  =  Ax  +  X(6,  x,  y,  z) 
y  =  By  +  Y(6,x,y,z) 
z  -  Cz  +  Z{d,x,y,z) 


(C.l) 


here  x,  y  and  z  are  vectors  in  for  n(x),  n(y)  and  n(z)  their  corresponding  dimen¬ 

sions. 

For  system  (C.l)  we  look  for  a  locally  invariant  manifold  parametrized  by  the 
variable  y,  such  a  manifold  is  known  as  the  center  manifold.  For  more  details  we  refer 
the  reader  to  [5]  and  [9]. 

We  make  the  following  assumptions; 

(hi)  A,B,C  are  constant  matrices  in  real  canonical  form.  A  has  eigenvalues  with 
negative  real  part,  B  has  eigenvalues  with  zero  real  part  and  C  has  eigenvalues  with 
positive  real  part, 

(h2)  X,y,  Z  are  vector  valued  functions  of  class  in  a  neighborhood 
Ns{{9,x,y,z)/  I  X  P  +  I  y  P  +  U  P< 

(h3)  X^Y^Z  and  their  derivatives  with  respect  to  $,  x,  y  and  z  are  all  zero  for 
arbitrary  values  of  9  when  (x,y,z)  =  0. 

(h4)  For  matrices  Rx,  Ry  and  Rx  of  the  form 


Ra  = 


with  Ii^a  and  72, «  identity  matrices  of  given  dimensions,  and  a  =  x,y,z,  we  have  the 
following  identities 

X{9  -b  l,x,y,z)  =  RxX{6,RxX,Ryy,Rxz) 

Y(6  +  l,x,y,2)  =  RyY{6,  RxX,Ryy,R,z) 

Z{9  +  1,  X,  y,  z)  =  RzZ(9,  R^x,  Ryy,  R~z) 
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and  RxA  =  ARx,  RyB  =  BRy  and  R^C  =  CRz, 

(h5)  For  some  constants  c  >  0,  p  >  0  and  >  0 

^(0,0, y,0)  <  c  I  y  |P,  ^(&,0,y,0)  <  c  ]  y  \p+\ 

|^(0,O,y,O)  <  c  I  y  |P,  ^(^,0,y,0)  <  c  |  y  |p+\ 

for  every  y  so  that  |  y  |< 

Remark  C.l  We  note  that  (h4)  implies  that  the  functions  X,  Y  and  Z  are  2-periodic 
in  d.  Then  we  can  apply  the  Center  Manifold  Theorem  to  system  (C.l).  We  want  to 
use  (h4)  and  (h5)  to  obtain  extra  properties  for  the  center  manifold. 

Theorem  C.l 

For  system  (C.l )  there  is  a  center  manifold  given  by 

M  =  {(0,x,y,z)/  0  arbitrary,  |  y  |<  ^2,  a:  =  u{9,y),  z  =  v(d,y)] 
where  the  functions  u  and  v  are  defined  and  of  class  in  a  neighborhood 

^sz  =  ((0,y)/  \y\<S2} 


and  they  satisfy 


u($,Ryy)  =  Rxti(0,y) 
v{0,Ryy)  =  Rzv{9,xj) 


{C.2) 


And  for  some  constants  K  and  S3 


1  I  </f  I  yr,  l|^(9,y)l</'f|yr’, 
I  gj(«,y)  I  <  A'  I  y  r,  I  |^(«,y)  l<  A I  y  T'. 


(C.3) 


for  all  I  y  |<  <^3. 

Proof.  This  theorem  is  proved  by  using  the  contraction  principle  in  an  adequate  Banach 
space.  See  [9].  We  start  describing  a  change  of  variables  designed  to  deal  with  the  local 
nature  of  the  problem.  Let  (f>{r)  be  a  smooth  real  valued  function  so  that  (f>{r)  =  1  if 
0  <  r  <  1/2,  and  =  0  if  1  <  r  <  00.  For  large  R  and  for  A  in  R  we  define 


Xie,  x,y,z,\)  =  <i>i\x\^  +  \y\'^  +R\'^)X~^  X{6,  Ax,  Ay,  \z) 

Y{e,x,y,z,\)  =  ^^(1  X  p  +  I  y  p  +  I  z  p  +ilA2)A-'F(0,Ax,Ay,A^) . 
Z(0,x,y,2,  A)  =  <^(1  X  p  +  I  y  P  +  I  2  p  +R\^)\-^Z{e,\x,\y,\z) 
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(C.4) 


With  this  change  of  variables  system  (C.l)  is  transformed  into 

e  =  i 

X  =  Ax  +  X(6,  x,y,  z.  A) 
y  =  By  +  Y{6,x,y,z,\) 
i  =  C2  +  Z{d^x,y,  z,  A) 

and  the  functions  X,  Y  and  Z  satisfy  the  following  properties: 

(h2)A  X,  Y  and  Z  are  defined  and  continuous,  and  for  each  A  fixed  they  are  of 
class  in  (0,  x,y,2), 

(hS)^  X,Y,Z  and  their  derivatives  with  respect  to  6,  x,  y  and  z  are  all  zero  for 
arbitrary  values  of  6  when  (a:,y,  2)  =  0.  And  ^  ^  ^^{X,Y,  Z)  — ^  0  uniformly  in 
(0, i,y, 2)  as  A  —>  0  for  0  <|  p  |<  2, 

(h4)A 

X(0  +  l,2:,y,2)  =  RxX{9,RiX,Ryy,R2z), 

Y{6  +  l,x,y,2)  =  RyY(6,RxX,Ryy,R^z), 

Z(d  +  1,  X,  y,  2)  =  R^Z(e,  RxX,  Ryy,  R^z), 

(hS);^  For  some  constants  c  >  0,  p  >  0  and  >  0  and  for  A  <  1 

|^(9,0,y,0.X)  <  c  I  y  r.  ^(«.0,y,0,A)  <  c  1  y  r-, 

^(S.o.y.o,  A)  <  c  I  y  p,  ^(«,0,y,0,A)  <  c  |  y  r‘. 

for  every  y  so  that  |  y  |<  <^1- 

We  define  the  following  space  of  functions 

=  {{u{6,xj),v{9,y))/u  and  v  satisfy(fl)  —  (f4)} 

(fl)  u  and  V  are  of  class  for  all  6  and  y, 

(f2)  u  and  v  have  period  2  in  6, 


(f3)  u(^,0),  v(e,0),  §f(9,0)  §^(9,0)  are  all  zero, 
(f4) 


(u,t;)  11=  max  sup  |  Ju,v)  1<  oo. 
0<|p|<l(e,y)  ^ 


R}  provided  with  the  norm  given  in  (f4)  is  a  Banach  space.  For  (u,  v)  S  R}  and  for 
^0)  yo  we  define  ^O)  yoi  X)  and  ^(t,  9q,  yo,  A)  as  the  unique  solution  of  the  differential 
equation 

^  ^  (C-5) 

y  =  By  +  Y{9,u{9,y),y,v{9,y),X) 
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and  for  which  ^(0,6/o,yo,  A)  =  6o  and  ‘n(0,9o,yo,X)  =  yo-  Obviously  the  function  V’  is 
given  simply  by  ip{t)  =  do+t-  We  define  now  our  operator  T  va.'H} .  The  u-component 
of  T  is  given  by  the  formula 

r° 

r„(u,u)(0o,yo)  =  /  e  ^‘"X{t()((T),u(xp((T),T){a)),T]((r),v{ip((r),7](cr)),X)da  ((7.6) 

J  —  OO 

and  similarly  the  i>- component  is  defined  by 

r° 

r„(u,t;)(^o,yo)  =  /  e  ^^Z{xp((r),u(tp(a),T]((T)),T]((r),v(ip(a),T](a)),X)da  (C.7) 

J  OO 

In  [9]  it  is  shown  that  for  A  small  enough  the  operator  T  is  a  contraction  in  the  unit 
ball  of  'HX  .  We  will  consider  the  following  closed  subset  of  : 

7^^(Jv )  =  {(u, u)  G  (^>^)  satisfies  (C.2)  and  (C.3)} 

In  what  follows  we  show  that  for  a  certain  constant  K  the  operator  T  is  a  contraction  in 
We  first  show  that  if  (u,u)  satisfies  (C.2)  then  (Tu(u,i;),T„(u,u))  also  satisfies 
(C.2).  We  claim  that 


rj{t,0Q  +  l,J/o,A)  =  Ryr](t,Oo,Rytjo,X)- 


((7.8) 


Lets  assume  for  the  moment  that  (C.8)  is  true.  Since  (u,u)  satisfies  (C.2),  by  (h4)A  and 
(C.8)  we  have 

Tv{u,v)(0q  +  l,yo)  = 


-f 


^  X{dQ  +  1  +  <7,  w(^o  +  I  +  <7, 7?(<T,  Oq  +  1,  yo)  A)), 


■n{a,9Q  +  l,yo,A)),j;(0o  +  1  +  cr,  ^(a,  0o  +  1,  yo,  A)),  A)  dcr  = 

fO 

=  e  '^’'RxX{9o  +  c,Rj;u{9o  +  1  +  cr,r}{a,9o + 

J  — OO 

Ryr]{a,9o  +  l,yQ,X)),Rzv{9o  +  l  +  a,r]{cr,  9o  +  l,yo,  A)),  A)  dcr 

/.O 

=  Rx  e  ^‘'X{9o  +  (r,u{9Q  +  <x,Ryr]{a,9o  +  l,yo,X)), 

J  —  OO 


r]{a,9o,Ryyo,X)),v{9Q  +  <T,Ryr]{a,9o  +  1,  yo,  A)),  A)  da 


=  Rr  f 
J  —  o 


e  ^"X{9o  +(7,u(9q  +  a,r](cx,9o,Ryyo,X)), 


T;((7,do,i2yyo,  A)),u(do  +  (T,T](cr,9o,Ryyo,X)),X)dcr 
=  RxT{u,v){9o,Ryyo) 


In  other  words  we  have  proved 


ru(«,t;)(^o  +  l,yo)  =  RxTu{u,v)(0Q,Ryyo)  (C.9) 

In  a  similar  feishion  we  can  show 

Tv{u,v)(6o  +  l,yo)  =  RzTv(u,v)(6o,Ryyo)  (C.IO) 

In  order  to  complete  the  proof  we  check  that  (C.8)  is  correct.  Let  us  consider  the 
function 

T]*(t)  =  Ryr)(t,eo,Ryyo). 

Then  77*(0)  =  yo,  and  since  (u,w)  satisfies  (C.2)  and  Y  satifies  (h4)A  we  have 

=RyBr](t,  6(i,RyyQ )  +  RyY(6Q  +  t,  u(9o,r}),  rj,  v(0q  +  t,  y),  A) 

=Bt]*  +  F(5o  +  1  +  Rxu{6q  + 1, 7?),  Ryt],  Rzv{6o  + 1 , »?),  A) 

=Brj*  +  Y{6q  +  1  +  i, u{9q  +  1  t,T]*'),y* ,v{6o  +  1  +  t)^*))A) 

Then,  by  uniqueness  of  the  solution  of  (C.5)  we  have 

T]*{t)  =  y{t,6Q  +  l,yo,  A), 


i.e.  (C.8)  is  satisfied. 

We  now  show  that  a  constant  K  exists  so  that  given  (u,  u)  G  that  satisfies  (C.3) 
the  functions  (r„(u,y),T„(u,u))  also  satisfies  (C.3).  We  will  use  the  following  result  on 
the  behaviour  of  ri{t,  0Q,yo,\)  that  can  be  prove  using  the  ideas  given  in  [9],  Lemma  3: 
for  A  small  enough  and  (u,v)  €  W  we  have 

U((.«o,!/o,A)|<e^’l'l|!/ol  (C.12) 


I  ^(i,9o.!/o,A)|<e"’W,  (C.13) 

I  |<  (1+  I  i  I  !/„  I  .  (C.14) 

The  constant  7  has  been  chosen  so  that  -a4-27(p+2)  <  0  where  a  and  —a  are  lower  and 
upper  bound  for  the  real  part  of  the  eigenvalues  of  the  matrices  A  and  C  respectively. 

We  have  to  find  K  so  that  if  (u,u)  satisfies  (C.3)  then  T{u^v)  also  satisfies  (C.3). 
Let  9{0Q,yo)  =  T„(«,u)(^o,yo),  then  we  have 


{C.16) 


From  (C.3)  and  (C.12)  we  have 


du 

dy 


(00  +a,7,((T,0o,yo,A))  1<  K  1 7,(a,0o,yo,A)  K<  1  yo  r, 


and  by  integrating  (C.3)  we  obtain 


u{do,yo)  I< 


K 

p+1 


VO  1'’+' 


and 


v(0o,yo) |< 


K 

P+1 


yo 


|p+i 


(ai7) 


(C.18) 


Let  us  estimate  the  middle  term  in  (C.15).  Noting  that  X  is  of  class  C^,  by  the  Mean 
Value  Theorem  and  using  (h3)A,  we  see  that  for  A  small  enough 

I  ^(^0  +  cr,  u,  y,  u,  A)  -  ^(^0  +  0,  y,  0,  A)  |< 

<  ^(^0  +^^,y(^,^o,yo,A))  I  +  1  u(0o  +iT,7?(a,0o,yo,A))  1}  (C.19) 

<  K  I  T}((T,9o,yQ,\)  1^+^ 

<  AV^(p+1)I‘^I  (  yo  |P+\ 

here  we  used  (C.12),  (C.17)  and  (C.18).  Then,  by  (h5)x,  (C.13)  and  (C.19)  we  obtain 

I  r  e(“-2^>"(c  +  A^e-2'^(^+')"|yo|)dt}|yor  (C.20) 

7-00  oy  7-00 

Choose  K  so  that  ^ 

f  <  E  (C.21) 

7-00  6 


and  choose  64  so  that. 


i-27(p+2))<t^^  < 

-  6 


ro 

J  —00 

Then,  from  (C.20),  (C.21)  and  (C.22)  we  conclude  that  for  all  |  yo  1<  84 

r«  _.„dXdrj 


(C.22) 


(C.23) 


Now  we  estimate  the  first  term  in  the  right  hand  side  of  (C.15).  By  (C.13)  and  (C.16) 


/ 


dx  dy  dy 


/O  PiY 

^(„-27(p+1))<./,'  I  _  I  I  |P  (C.24) 


From  (h3)A  and  by  restricting  A  if  necessary,  we  have 


r  Ja-2t(p+l))a  I  ^ 
7-00  '  dx 


da  <\ 
~  3 


and  then  from  (C.24)  and  (C.25) 


^-A<T 

dx  dy  dy 


Uo  r 


In  an  analogous  fashion  we  can  prove  that 


dz  dy  dy 


yo  r . 


(C.25) 


(C.26) 


(C.27) 


Thus,  from  (C.18),  (C.23),  (C.26)  and  (C.27)  we  see  that  there  is  <^5  >  0  and  K  so  that 
for  all  I  j/n  |<  ^5 

||^(^o,yo)|</i  |yol".  (C.28) 

A  similar  argument  can  be  given  to  show  that  if  h[6o,yo)  =  T'„(u,u)(0o,  J/o)  then 


|^(^o,yo)|<A'lyor.  (C.29) 

In  order  to  show  that  the  derivative  with  respect  to  0  of  ^  and  h  satisfies  condition  (C.3) 
we  proceed  in  a  similar  fashion.  We  may  need  to  reduce  further  A  and  increase  K.  We 
omit  the  details.  0 
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